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PREFACE 


This  report  describes  studies  of  the  phase  composition  of  pore  water  in  cold  rocks.  It 
forms  part  of  a  broader  investigation  covered  by  DA  Task  1T062112A13001,  Cold  Regions 
Research-Applied  Research  and  Engineering,  Work  Unit  Mechanics  of  Frozen  Rocks.  The 
work  was  directed  by  Dr.  Malcolm  Mellor,  Research  Civil  Engineer  of  the  Applied  Research 
Branch  (Mr.  A.F.  Wuori,  Chief),  Experimental  Engineering  Division  (Mr.  K.A.  Linell,  Chief). 
The  report  has  been  reviewed  by  Dr.  Pieter  Hoekstra,  Physical  Sciences  Branch,  Research 
Division,  USA  CRREL. 

The  writer  gratefully  acknowledges  the  help  and  guidance  of  Dr.  Duwayne  Anderson 
and  Dr.  Pieter  Hoekstra,  and  the  active  assistance  and  experimental  contributions  of 
Mr.  Alan  Greatorex  and  Specialists  Russell  Rainey,  Ronald  Severson,  and  Richard  Testa. 

The  contents  of  fliis  report  are  not  to  be  used  for  advertising,  publication,  or  promotional 
purposes.  Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval 
of  the  use  of  such  commercial  products. 


Ill 


CONTENTS 

Page 

Preface .  ii 

General  introduction .  1 

A.  Pore  water  freezing  data  from  differential  thermal  analysis .  2 

Apparatus  and  technique .  2 

Results .  3 

Discussion . 7 

Conclusions .  14 

B.  Indirect  determination  of  pore  water  freezing  data  for  rocks .  14 

Air  penetration  tests . 15 

Mercury  penetration  r$e.«surements  .  15 

Calculation  of  unfrozen  water  content  and  freezing  point  depression  17 

Comparison  of  calculated  and  measured  freezing  characteristics  .  19 

Conclusions.. .  21 

C.  Electrical  conductivity  measurements  .  21 

Procedures . 22 

Results .  24 

Discussion  of  results .  24 

D.  Thermal  strains  in  cold  rock  . 31 

Preliminary  tests  . 31 

Recording  dilatometer .  32 

Test  procedure  .  36 

Results .  37 

Discussion .  40 

Conclusions .  43 

E.  Isothermal  compressibility  of  cold  rocks .  43 

Test  method .  44 

Test  results .  46 

Discussion  of  results  .  46 

Conclusion  and  recommendations .  50 

General  summary  of  results  .  50 

Literature  cited .  51 

Appendix  A:  Water  adsorption  .  55 

Appendix  B:  Adsorption  and  absorption  by  Rochester  shale  .  57 

Appendix  C:  Low  temperature  conductivity  of  saturated  wood  .  59 

Abstract .  61 

TABLES 

Figure 

1.  Details  of  DTA  cell .  3 

2.  Experimental  set-up .  4 

3.  Freezing  point  depression  as  a  function  of  water  content  for  Berea 

sandstone .  8 

4.  Freezing  point  depression  as  a  function  of  water  oontent  for  Indiana 

limestone .  8 


IV 


r  ' 


CONTENTS  (Cont’d) 

ILLUSTRATIONS  (Cont’d) 

1  Page 

5.  Logarithmic  plot  of  freezing  point  depression  against  water  content 

for  Berea  sandstone  . 9 

li.  Logarithmic  plot  of  freezing  point  depression  against  water  content 

for  Indiana  limestone .  9 

7.  Supercooling  as  a  function  of  water  content  for  Berea  sandstone .  10 

cS.  Supercooling  as  a  function  of  water  content  for  Indiana  limestone .  10 

9.  Latent  heat  as  a  function  of  water  content  for  freezing  of  Berea  sandstone  11 

10.  Latent  heat  as  a  function  of  water  content  ft  r  freezing  of  Indiana  limestone  1 1 

11.  Logarithmic  plot  of  latent  heat  against  water  content  for  freezing  of  Berea 

sandstone .  12 

12.  Logarithmic  plot  of  latent  heat  against  water  content  for  freezing  of 

Indiana  limestone .  12 

10.  Air  penetration  data  for  Berea  sandstone  and  Indiana  limestone  .  lb 

11.  Air  penetration  data  for  Barre  granite .  lb 

lh.  Mercury  penetration  data  for  Berea  sandstone .  18 

lb.  Mercury  penetration  data  for  Barre  granite  .  18 

17.  Freezing  point  depression  as  a  function  of  water  content  for  Berea 

sandstone .  19 

18.  Freezing  point  depression  as  a  function  of  water  content  for  Indiana 

limestone .  19 

19.  Freezing  point  depression  as  a  function  of  water  content  for  Barre  granite  20 

20.  Test  cell  used  for  low  temperature  electrical  conductivity  measurements  28 

ill.  Arrangement  of  test  cell  in  Dewar  flask  .  23 

22.  Conductivity  at  104  Hz  as  a  function  of  temperature  for  Barre  granite  . .  25 

23.  Conductivity  at  104  Hz  as  a  function  of  temperature  for  Indiana  limestone  25 

21.  Conductivity  at  104  Hz  as  a  function  of  temperature  for  Berea  sandstone  2b 

25.  Conductivity  at  104  Hz  as  a  function  of  temperature  for  fine-grained 

polycrystalline  He .  26 

2b.  Conductivity  at  104  Hz  as  a  function  of  temperature  for  fine-grained 

polycrystalline  ice .  27 

27.  Conductivity  plotted  against  1/T  for  Barre  granite  .  27 

28.  Conductivity  plotted  against  1/T  for  Indiana  limestone .  28 

29.  Conductivity  plotted  against  1/T  for  Berea  sandstone  .  28 

30.  Conductivity  plotted  against  1/T  for  fine-grained  polycrystalline  ice  . . .  29 

31.  Optical  ext ensometer  results  for  air-dry  Berea  sandstone  .  32 

32.  Optical  extensometer  results  for  air-dry  Indiana  limestone .  32 

33.  Optical  extensometer  results  for  air-dry  Barre  granite .  32 

31.  Recording  dilatometer  .  ' ' 

35.  Recording  dil at ometer  . 

3b.  Calibration  data  from  recording  dilatometer .  *■> 

37.  Calibration  data  from  optical  extensometer  .  '*•> 

38.  Calibration  data  from  optical  extensometei  .  3b 

29.  Thermal  strain  iceoids  for  Berea  sandstone  . 

Hi.  Thermal  strain  lecmds  loi  Indiana  limestone .  ’  ^ 


1 


V 


CONTENTS  (Cont’d) 

ILLUSTRATIONS  (Cont’d) 

Figure  Page 

41.  Thermal  strain  records  for  Barre  granite  .  :{8 

42.  Coefficient  of  linear  expansion  plotted  against  water  content  for  Berea 

sandstone  and  Indiana  limestone .  j:i 

43.  Phase  diagram  for  ice  in  the  pressure-temperature  plane .  U 

44.  High  pressure  die  and  loading  machine  .  15 

45.  Pressure-volume  relationships  for  hydrostatic  compression  of  dry  sand¬ 

stone  and  water-saturated  sandstone  at  -10  C  .  47 

46.  Volumetric  strain  as  a  function  of  hydrostatic  pressure  at  -10  C  for  dry 

sandstone  and  water-saturated  sandstone  .  18 

47.  Pressure-volume  relationships  for  hydrostatic  compression  of  nominally 

dry  and  nominally  saturated  specimens  of  Indiana  limestone  at  -103C  49 

48.  Pressure-volume  relationships  for  hydrostatic  compression  of  nominally 

dry  and  nominally  saturated  specimens  of  Barre  granite  at  -10°C  ....  49 

TABLES 

Table 

I.  DTA  results  for  Berea  sandstone .  4 

II.  DTA  results  for  Indiana  limestone .  5 

III.  DTA  results  for  Berea  sandstone .  6 

IV.  DTA  results  for  Barre  gray  granite  .  7 

V.  General  properties  of  test  specimens  .  22 

VI.  Thermal  strain  data  for  cold  rocks  .  39 

VII.  Linear  freezing  strains .  42 

VIII.  Compressibility  of  ice  at  -10°C .  50 


BLANK  PAGE 


PHASE  COMPOSITION  OF  PORE  WATER  IN  COLD  ROCKS 

by 

Malcolm  Mellor 


GENERAL  INTRODUCTION 

In  the  past,  little  systematic  attention  has  been  given  to  the  effects  of  cold  environments  on 
rocks  (as  distinct  from  soils),  and  thus  a  study  was  undertaken  at  USA  CRREL  to  determine  how 
mechanical  properties  of  rocks  are  affected  by  cold  environments,  and  in  particular  by  freezing  of 
interstitial  water.  Since  it  was  anticipated  that  mechanical  properties  of  cold  rocks  would  be 
influenced  significantly  by  changes  in  state  of  the  pore  water,  a  subsidiary  investigation  into 
phase  changes  of  pore  water  at  temperatures  below  0°C  was  also  undertaken. ;  That  subsidiary 
investigation  is  the  subject  of  this  report. 

The  main  expectation,  derived  largely  from  the  literature  on  properties  of  frozen  soils,  was 
that  freezing  of  pore  water  would  be  controlled  by  pore  geometry  and  interface  properties,  so  that 
some  water  would  remain  liquid  at  temperatures  well  below  0°C.  Another  expectation,  prompted 
by  the  phase  diagram  for  bulk  ice,  was  that  frozen  pore  water  might  undergo  phase  transformation 
at  constant  temperature  if  sufficiently  high  pressures  were  reached. 

The  general  problem  was  a  twofold  one:  1)  to  determine  the  state  of  pore  water  as  a  function 
of  temperature  at  constant  pressure,  and  2)  to  determine  the  state  of  pore  water  as  a  function  of 
pressure  at  constant  temperature. ;  The  first  requirement  was  for  a  determination  of  the  freezing 
characteristics  of  the  rocks  used  in  the  mechanical  test  program,  i.e.  it  was  necessary  to  estab¬ 
lish  the  temperature  for  initial  nucleation  of  ice  and  the  relationships  between  liquid  water  content 
and  temperature. ;  Other  matters  of  interest  were  possible  discontinuities  in  the  temperature  char¬ 
acteristics  for  pore  water,  and  direct  mechanical  effects  connected  with  thermal  strain.  The 
investigation  of  pressure  effects  was  concerned  mainly  with  the  possibility  of  transformations  from 
Ice-I  to  water  and  to  high  pressure  modifications  of  ice  in  the  frozen  pore  water. 

Since  the  study  had  to  run  concurrently  with  the  main  program  of  mechanical  tests,  ready 
availability  of  techniques  and  equipment  was  an  important  practical  consideration. 

Several  distinct  experimental  techniques  were  employed,  and  the  following  report  consists  of 
five  separate  parts  describing  these  techniques  and  the  results  obtained  from  them.  The  subjects 
covered  are: 

A.  :  A  determination  of  pore  water  freezing  data  by  differential  thermal  analysis. 

B.  A  study  of  indirect  (room  temperature)  determination  of  pore  water  freezing  characteristics. 

C.  Measurements  of  electrical  conductivity  as  a  function  of  temperature. 

D.  Measurements  of  the.'mai  strain  at  low  temperatures. 

E.  Isothermal  compressibility  measurements  at  -10°C. 
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'Hirsc  studies  were  undertaken  as  part  of  an  engineering  investigation,  and  not  as  fundamental 
research .  Nevertheless,  they  appear  to  lie  highly  relevant  to  geological  problems  of  frost  weathering 
and  to  problems  of  frost  deterioration  in  synthetic  porous  materials  such  as  concrete,  bricks  and 
ceramics. 


A.  FORE  WATER  FREEZING  DATA  FROM  DIFFERENTIAL  THERMAL  ANALYSIS 

In  stinking  phase  composition  of  pore  water  in  rocks  subjected  to  low  temperatures,  it  was 
first  necessary  to  determine  the  temperature  at  which  freezing  of  pore  water  is  initiated  for  a  given 
water  content,  and  to  determine  the  amount  of  water  which  actually  changes  to  the  solid  state  at 
this  temperature.  The  method  used  was  differential  thermal  analysis  (DTA),  in  which  the  test 
sample  and  an  inert  reference  material  are  heated  or  cooled  at  the  same  rate  while  the  temperature 
difference  between  sample  and  reference  is  recorded.  A  phase  change  in  a  component  of  the  test 
sample  produces  a  perturbation  in  the  record  of  temperature  difference  versus  time.  The  rocks 
studied  were  Barre  granite,  Indiana  limestone,  and  Berea  sandstone;  their  water  contents  were 
varied  from  zero  to  complete  saturation. 

Apparatus  and  technique 

Differential  thermal  analysis  lias  been  used  to  study  the  freezing  of  pore  water  in  rocks  by 
Dunn  and  Hudec  (1965a,  b),  wtio  built  a  precise  apparatus  and  calibrated  it  to  give  absolute  values 
of  heat  output  during  the  freezing  of  pore  water.  For  the  present  work  a  simpler  arrangement  was 
used;  no  attempt  was  made  to  obtain  an  absolute  calibration  for  output  of  latent  beat,  since  the 
method  seems  inherently  unsuitable  for  precise  quantitative  measurements. 

The  system  used  was  basically  similar  to  one  employed  at  the  U.S.  Army  Cold  Regions  Re¬ 
search  iuid  Engineering  Laboratory  by  Anderson  and  Tice  (personal  communication)  foi  studies  on 
frozen  soils.  Sample  and  reference  are  housed  in  a  metal  cylinder,  which  is  suspended  above  a 
pool  of  liquid  nitrogen  in  a  Dewar  flask  and  cooled  by  nitrogen  vapor  circulated  by  a  constant  speed 
fan.  Two  thermocouples  are  embedded  in  the  test  sample,  one  linked  to  another  thermocouple  in 
the  reference  material  to  give  differential  temperature  readings,  and  a  second  one  linked  to  an  ice- 
bath  junction  to  give  the  actual  temperature  of  the  sample.  Outputs  from  the  two  sensors  are  record¬ 
ed  simultaneously  on  a  two-channel  potentiometrie  recorder.  The  signal  from  the  differential  tempera¬ 
ture  sensor  is  passed  through  an  amplifier  before  reaching  the  recorder. 

The  first  DTA  cell  to  be  built  was  simply  a  scaled-up  version  of  the  cell  used  bv  Anderson  and 
Tice  for  soil  studies.  It  was  a  2-in. -diameter  cylinder  of  low  conductivity  stainless  steel,  23A  in. 
long,  bored  out  at  one  end  to  accommodate  a  rock  sample  1  in.  in  diameter  by  1  in.  long.  The  sample 
end  was  covered  by  a  H-in. -thick  stainless  steel  cap,  and  the  reference  material  was  the  solid  end 
of  the  main  cylinder.  This  cell  was  unsatisfactory  because  of  poor  matching  of  thermal  conductivity 
between  the  stainless  steel  reference  and  the  rock.  When  modified  with  a  Lucite  plug  at  the  refer¬ 
ence  end,  there  remained  objections  on  grounds  of  lack  of  symmetry. 

The  final  design  of  the  cell  was  that  shown  in  Figure  l.  The  test  sample  and  the  reference 
plug,  each  1  in.  diameter  by  1  in.  long,  were  housed  at  opposite  ends  of  a  2-in. -diameter  copper 
cylinder,  which  was  closed  at  both  ends  by  \-in.-thiek  copper  cover  plates.  The  copper  cylinder, 
the  rock  samples,  and  the  glass  reference  plug  were  all  machined  precisely  to  fit  within  0.002  in. 

'I  lie  thickness  and  the  thermal  conductivity  of  the  copper  cylinder  were  sufficient  to  damp  out 
spatial  and  temporal  fluctuations  o(  surface  temperature  resulting  from  imperfect  circulation  of  the 
cool  mg  va|wii . 
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Thermocouples  were  placed  at  the  center  of  the  test  sample 
and  the  reference  plug.  A  hole  in.  in  diameter  was  sufficient 
to  admit  two  pairs  of  double-insulated  T3  gauge  thermocouple 
wires.  The  thermocouples  were  cemented  in  place  with  cold- 
set  lint;  epoxy  adhesive,  which  was  carefully  worked  m  to  eiisute 
that  all  remainin';  space  in  the  thermocouple  hole  was  tilled.  The 
thermocouples  were  copper-const  ant  an,  and  the  actual  junctions 
were  trimmed  to  a  length  of  in.  The  holes  in  the  cover  plates 
throimh  winch  the  thermocouple  leads  passed  were  sealed  with 
silicone  rubber  adhesive. 

After  loading  at  room  temperature,  the  l)TA  cell  was  suspended 
over  liquid  nitrogen  in  the  Dewar  flask.  A  stirring  propeller  was 
positioned  near  the  cell,  and  the  flask  was  closed  with  a  g'-in.- 
t hick  split  cover  of  foamed  polyurethane.  Cooling  rate  was  adjust¬ 
ed  roughly  by  proximity  of  the  cell  to  the  liquid  nitrogen  surface 
and  by  the  stirring  rate:  the  cell  was  typically  about  1  m.  above 
the  liquid  nitrogen,  and  stirring  rate  was  normally  maintained  at 
700  rpm.  Cooling  rates  between  0  and  -10  C  were  mainlv  in  the 
range  1  -  2  C  min. 

In  principle  the  apparatus  should  have  been  capable  of 
detecting  temperature  differences  between  sample  and  reference  of 
0.001  C,  but  the  electronic  noise  level  was  such  that  only  tempera¬ 
ture  differences  in  excess  of  0.01  C  could  he  detected.  This 
performance  could  be  improved  significantly  with  further  work  on 
the  electronic  equipment. 

Each  rock  type  rested  was  loaded  into  the  cell  initially  while  still  warm  after  oven  drying. 

The  sample  was  run  dry  to  check  the  apparatus  and  the  scale  settings,  and  to  establish  a  base  line. 
Distilled  water  was  added  to  the  surface  of  the  high-porosity  samples  from  a  vernier  syringe.  The 
water  soaked  into  the  rock  quite  readily,  but  it  apparently  took  many  hours  to  diffuse  throughout 
the  sample.  A  cycle  of  cooling  and  warming  through  10  C  or  more  speeded  up  the  diffusion  of 
newly  introduced  water.  To  obtain  complete  saturation  the  entire  cell  was  placed  in  a  vacuum 
vessel,  and  the  sample  was  removed  sufficiently  to  permit  immersion  in  distilled  water  without  dis¬ 
connecting  the  wiring. 

There  is  a  possibility  that  some  water  was  adsorbed  from  the  air  before  the  first  increment  of 
water  from  the  syringe  was  taken  up. 

Figure  2  shows  the  experimental  set-up,  which  includes,  from  left  to  right,  a  stabilized  d-e 
microvolt  amplifier  (Leeds  and  Northrup  98135)  and  a  potentiometric  recorder  (Moseley  710013),  an 
ice  bath,  the  main  Dewar  flask  and  stirring  motor,  and  a  constant  speed  control  unit  (Cole-Partner 
model  44130). 

Results 

Basic  results,  in  the  form  of  clout  records,  are  not  reproduced  here,  but  data  abstracted  tiem 
the  charts  are  listed  in  Tables  I-IV.  Values  given  for  supercooling  are  the  temperatures  at  which 
nuclealioti  occurred.  Values  given  for  freezing  point  depression  are  the  temperatures  to  which  the 
samples  rose  after  initial  nucleation.  Values  given  for  latent  heat  (in  arbitrary  units)  are  obtained 
from  the  areas  of  the  perturbations  of  temperature  difference,  i.e.  they  give  a  measure  of  the 
integral  of  temperature  difference  with  respect  to  lime  for  the  duration  of  the  freezing  pulse. 
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Figure  I.  Details  of  DTA 
cell. 


4 


PHASE  COMPOSITION  OF  P&BE  WATER  IN  OOLD  ROCKS 


Figure  2.  Experiment*/  set-up.  The  items  shown  are,  from  left  to  right, 
an  amplifier  surmounted  by  a  potentiometric  recorder,  an  ice  bath,  a  Dewar 
flask  containing  the  cell  and  liquid  nitrogen,  with  fan  motor  above  and  a 
constant  speed  control  unit  for  the  fan. 


Table  I.  DTA  racilta  for  Bant  aacdctowe  (drat  tariaa  of  tecta). 

Sample  dimensions:  1.010  In.  in  diameter  x  0.998  in.  in  length.  Axial  thermocouple  hole  0.003  in.  in 
diameter  x  0.52  in.  in  length. 

Sample  volume  (adjusted  for  thermocouple  bole):  0.797  in.1  ■=  13.08  cm* 

Sample  mass:  27.7  g 

Dry  bulk  density  of  sample:  2. 1 19  g/cm* 

Effective  porosity:  19.87. 

Vaier  added  to  sample:  Distilled  deionised  water,  conductivity  10**  (ohm-cm)*'.at  26°C  and  60  Hz. 


Water 

Vater  content 

content  g  water/ 

Test  g  g  rock 


Pressing 

Super-  point  Latent  heat  Cooling 

codling  depression  arbitrary  rate 

°C  °C  units  °C/min 


1 

0 

0 

2 

0 

0 

3 

0.2 

7.22  x  10*’ 

2.76 

1. 14 

0.40 

1.17 

4 

0.2 

7.22  x  10*' 

3.33 

1.61 

0.40 

1.66 

5 

0.4 

1.44  x  10’* 

2.15 

0.43 

0.63 

1.14 

6 

0.4 

1.44  x  10*’ 

2.60 

0.41 

0.58 

0.89 

7 

0.6 

2.16  x  10’1 

2.08 

0.35 

0.60 

1.40 

8 

0.6 

2.16  x  10"* 

2.49 

0.30 

0.61 

1.12 

9 

1.0 

3.61  x  10'1 

2.08 

0.26 

0.94 

1.04 

10 

1.0 

3.61  x  10*’ 

2.47 

0.30 

0.84 

2.00 

11 

1.5 

5.41  x  10*1 

1.78 

0.26 

1.50 

1.16 

12 

1.6 

5.41  x  10*1 

2.34 

0.27 

1.28 

1.49 

13 

2.4 

8.66  xHD*1 

0.17 

0.08 

2.06 

1.36 

14 

2.4 

8.66  x  10'* 

1.42 

0.16 

2.03 

1.68 

15 

2.4 

8.66  x  10'* 

0.88 

0. 12 

2.13 

1.43 
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Table  D.  DTA  results  for  Indiana  (Bedford)  limestone. 

Sample  dimensions:  1.010  in.  in  diameter  x  0.998  in.  in  length.  Axial  thermocouple  hole  0.063  in.  in 
diameter  x  0.52  in.  in  length. 

Sample  volume  (adjusted  for  thermocouple  hole):  0.797  in.1  13.08  cm1 

Sample  mass:  29.95  g 

Dry  bulk  density  ot  sample:  2.292  g/cm1 

Effective  porosity:  14.1% 

W'ater  added  to  sample:  Distilled  deionized  water,  conductivity  10"*  (ohm-cm)’1  at  25°C  and  60  Hz. 


Water 

Freezing 

Water 

content 

Super- 

point 

Latent  heat 

Cooling 

content 

g  water/ 

cooling 

depression 

arbitrary 

rate 

Test 

e 

g  rock 

°C 

°C 

units 

°C/min 

26 

0.01 

3.34  x  10'4 

27 

0.01 

3.34  x  10’4 

28 

0.02 

6.68  x  10'4 

29 

0.02 

0.68  x  10'4 

30* 

0.04 

1.34  x  10'5 

(9.01) 

(9.01) 

(0.03) 

1.19 

31 

0.04 

1.34  x  10's 

32* 

0.06 

2.00  x  10‘s 

(8.77) 

(8.77) 

(0.03) 

1.55 

33 

0.06 

2.00  x  10‘J 

34 

0.08 

2.07  x  10'1 

8.42 

8.42 

0.052 

1.48 

35 

0.08 

2.67  x  10"1 

8.53 

8.53 

0.039 

1.92 

36 

0.08 

2.07  x  10- 

8.66 

8.66 

0.048 

1.90 

37 

0.10 

3.34  x  10- 

8.26 

8.26 

0.076 

2.82 

38 

0. 10 

3.34  x  10- 

8. 10 

8.16 

0.092 

2.71 

39 

0. 10 

3.34  x  10- 

8.74 

8.74 

0.090 

2. 19 

40 

0.12 

4.00  x  10- 

8.40 

8.05 

0.12 

1.40 

41 

0. 12 

4.00  x  10- 

8.79 

8.77 

0.12 

2.35 

42 

0. 12 

4.00  x  10- 

8.66 

8.53 

0.11 

1.58 

43 

0.25 

8.34  x  10- 

6.41 

4.22 

0.45 

1.47 

44 

0.25 

8.34  x  10- 

6.65 

4.28 

0.49 

2.21 

45 

0.25 

8.34  x  10- 

6.67 

3.80 

0.46 

1.67 

46 

0.50 

1.67  x  10- 

5.36 

1.58 

pegged 
(0.76  +  A) 

1.39 

47 

0.50 

1.67  x  10- 

5.88 

2.15 

0.92 

2.22 

48 

0.50 

1.67  x  10- 

7.02 

2.70 

0.86 

1.22 

49 

1.00 

3.34  x  10- 

6.54 

0.85 

pegged 

(1.61  +  A) 

1.03 

50 

1.00 

3.34  x  10- 

2.60 

1.08 

1.75 

1.23 

51 

1.00 

3.34  x  10- 

5.75 

0.95 

1.66 

1.70 

52 

1.0 

5.34  x  10- 

0.93 

0.51 

2.69 

1.19 

53 

1.6 

5.34  x  10- 

1.27 

0.77 

2.73 

1.97 

54 

1.6 

5.34  x  10- 

1.03 

0.59 

2.84 

1.46 

•Premature 
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Table  Ul.  DTA  results  for  Berea  sandstone  (second  series  of  tests). 

Sample  dtmensiii/is  l.ulo  111.  111  diameter  ■  0,1(08  in.  in  lens'tli.  Axial  llieruiornuple  hoi o  0.0(18  in.  in 
diameter  ■  0.58  in.  in  li'iifOli 

Sample  \  ol iirjif  (.idpiMed  fur  tliernincnuple  hide):  0.707  in.1  18.08  cm5 
Sample  mass  87. p 
D n  hulk  den  sire  of  sample  8,  110  n  cm’ 


fif/ei  me  ptirn.sirv  10. S' 


H  .107 

added  to  sample 

Di sti  1  led  deionized 

uatei,  conductivity  10*''  (olim-em)* 

al  85' C  and  00  Hz. 

Hater 

Freezing 

Haler 

eonleiil 

Super 

point 

Ln ten  1  heat 

Cooling 

content 

g  Will (7 

<<)<)!  1  lift 

depression 

arbitrary 

rate 

Test 

8’ 

A’  rock 

‘  C 

"C 

twits 

"C  / min 

111! 

0.01 

8.01  ■ 

10  4 

(11 

o.or> 

l.SOfi 

■  Id"’ 

ill.' 

o.or. 

1.805 

v  10  ’ 

(18* 

0. 10 

8.(11  • 

10  1 

(8.52) 

(0.085) 

(VI 

0.  10 

8.01  - 

1(1** 

i>r> 

0.  10 

8.01  « 

10"’ 

(Hi 

0.  Ui 

5.45  ■ 

10*’ 

8. GO 

8.44 

0.044 

1.01 

0? 

o.  u> 

5.45  • 

10*’ 

5.09 

0.086 

8.08 

(IS 

0.  U: 

5.45  • 

I0M 

4.88 

0.058 

1.55 

ISO 

0.  Ui 

5.45  > 

10"’ 

4.17 

0.058 

8.24 

70 

0.80 

7.88  ' 

10  1 

O  »JQ 

*)  •)') 

KilS'I 

0.18 

1.00 

;i 

0.80 

7.88  ■ 

10'* 

8.47 

2.47 

0. 18 

1.24 

7:.’l 

0. 10 

8.01  ‘ 

10"’ 

78i 

0,80 

7.88  ■ 

10"’ 

8.89 

1.01 

0.18 

1.10 

74 ' 

0.80 

7.88  * 

10"’ 

o  r.*» 

8,|lls 

8. 18 

0.17 

0.91 

7M 

0.80 

7.88  • 

10"’ 

8.00 

2.28 

0. 10 

1.79 

70 

0.10 

1.44  ■ 

10 

8.  10 

1.48 

0.80 

1.14 

77 

0,10 

1.41  ' 

1(1  ’ 

8.41 

1.S1 

0.88 

1.10 

78 

0.  to 

1.11  ' 

.1 

10 

8.00 

1.10 

0.81 

1.89 

70 

0.80 

8.8H  • 

.1 

10 

1.08 

0.88 

0.07 

1.48 

80 

0.80 

8.SH  ■ 

.5 

10 

8.88 

0.40 

0.75 

1.88 

81 

0.80 

8. 88  ■ 

.1 

10 

8.84 

0.88 

0.00 

1.08 

88 

1.80 

4.88  ■ 

.! 

10 

1.08 

0.40 

0.98 

8.77 

88 

1  80 

1.88  • 

10 

1.40 

0.85 

1.08 

1.89 

84 

1.80 

4.88  v 

10'* 

1.70 

0.74 

1.10 

1.40 

80 

1.00 

5.77  • 

10*’ 

8.84 

0.85 

1.89 

1.18 

8(1 

1.00 

5.77  ■ 

.! 

10 

8.47 

0.10 

1.45 

0.89 

87 

1.00 

.5.  i  i  • 

10"’ 

8.48 

009 

1.41 

1.08 

88 

8  .  80 

0.85 

.1 

10 

0.80 

0.80 

8.80 

1.85 

SO 

\.Yl 

0.85 

10' ' 

0.85 

0.85 

8.48 

1.50 

00 

..o 

0.85 

1(1 J 

0.17 

0.10 

8.88 

1  11 

'  I'll 

a  Ini  i 

('In  i  I 
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Table  IV.  DTA  results  for  Bane  gray  granite. 

Sample  dimensions :  1.010  in.  in  diameter  x  0.999  in.  in  length.  Axial  thermocouple  hole  0.003  in.  in 
diameter  x  0.51  in.  in  length. 

Sample  volume  (adjusted  for  thermocouple  hole):  0.798  in.5  13.09  cm5 

Sample  mass:  34.5  g 

Dry  bulk  density  of  sample:  2.636  g/cm5 

Effective  porosity:  0.69% 

Water  added  to  sample:  Distilled  deionized  water,  conductivity  It)"4  (olim-cm)'1  at  25  C  and  00  Ilz. 


Water  Freezing 


Test 

Water 

content 

e 

content 

g  water/ 

g  ruck 

Super¬ 

cooling 

°C 

point 

depression 

°C 

Latent  heat 

arbitrary 

units 

Cooling 

rate 

C  min 

55 

0 

0 

56 

0.090 

2.61  x  10'5 

1.81 

1.55 

0. 15 

2.27 

57 

0.090 

2.61  x  10'5 

4.17 

2.73 

0. 14 

1.28 

58 

0.090 

2.61  x  10‘5 

3.88 

2.42 

0.15 

1.27 

59 

0.090 

2.61  x  10‘5 

3.83 

2.21 

0.16 

1.45 

A  number  of  warming  runs  were  made;  data  from  these  are  not  tabulated,  but  certain  character¬ 
istics  of  these  curves  and  the  curves  resulting  from  pilot  tests  were  studied. 

Data  for  Berea  sandstone  in  Table  I  are  not  reliable  for  the  lower  water  contents,  and  the  data 
in  Table  III  should  be  used  in  preference  to  these  first  results.  The  trouble  in  the  first  set  of  re¬ 
sults  apparently  stemmed  from  uncontrolled  adsorption  during  addition  of  the  first  increments  of 
water,  and  possibly  from  segregation  if  water  at  the  thermocouple,  which  may  have  been  sealed 
imperfectly.  Actually,  it  may  be  that  all  the  lower  values  of  water  content  given  for  Berea  sandstone 
and  Indiana  limestone  are  too  small,  since  there  could  have  been  adsorption  from  the  air  when  the 
cell  was  loaded  initially  and  when  it  was  opened  for  the  first  incremental  addition  of  water. 

Figures  3  and  4  give  freezing  point  depression  as  a  function  of  water  content  for  Berea  sand¬ 
stone  and  Indiana  limestone.  In  Figures  5  and  6  the  same  data  are  plotted  on  logarithmic  scales  in 
order  to  separate  the  points  obtained  for  low  water  contents. 

Figures  7  and  8  show  the  degree  of  supercooling  against  water  content  for  the  sandstone  and 
the  limestone. 

Figures  9  and  10  give  the  latent  heat  evolved  during  the  first  phase  change,  in  arbitrary  units, 
as  a  function  of  water  content.  In  Figures  11  and  12  the  same  data  are  plotted  on  logarithmic 
scales. 

Discussion 

The  apparatus  served  its  limited  purpose  satisfactorily,  but  in  its  present  form  it  is  unsuitable 
for  work  on  low  porosity  rocks,  since  water  content  of  the  rock  cannot  be  controlled  or  measured 
with  sufficient  precision.  Permanent  attachment  of  the  thermocouples  prohibits  weighing  as  a  means 
of  measuring  water  conti  ...  The  alternative  of  sensing  temperature  in  t  he  jacket  surrounding  the 
sample,  as  adopted  by  Dunn  and  Hudec,  also  has  certain  drawbacks  due  to  temperature  lag  and 
changes  of  thermal  conductivity  in  the  sample. 
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F igure  3.  Freezing  point  depression  as  a  function  of  water  content  for  Berea  sandstone. 


Figure  4.  Freezing  point  depression  as  a  function  of  temperature  for  Indiana  limestone. 
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Figure  5.  Logarithmic  plot  ol  freezing  point  depression 
against  water  content  for  Berea  sandstone. 


Figure  6.  Logarithmic  plot  of  freezing  point  depression 
against  water  content  for  Indiana  limestone. 
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Figure  7.  Supercooling  us  a  function  of  water  content  for  Bereu  sandstone. 
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Figure  S.  Supercooling  as  a  function  of  water  content  for  Indiana  limestone. 
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Figure  9.  Latent  heat  as  a  function  of  water  content  for  freezing  of  Berea  sandstone. 


Figure  10.  Latent  heat  as  a  function  of  water  content  for  freezing  of  Indiana  limestone. 
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Figure  11.  Logarithmic  plot  ol  latent  heat  against  water  content  for 
freezing  of  Berea  sandstone. 


Figure  12.  Logarithmic  plot  of  latent  heat  against  water  content  for 
freezing  ol  Indiana  limestone. 
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In  a  discussion  of  work  by  Lange  and  McKim  (1963),  Miller  (1963)  criticized  calorimetric 
determinations  of  freezing  point  depression  for  ceramic  blocks  on  the  grounds  that  latent  heat 
liberated  on  freezing  of  small  water  contents  would  be  insufficient  to  raise  the  temperature  of 
the  block  from  the  undercooling  temperature  to  the  equilibrium  freezing  temperature.  Neglecting 
the  effects  of  conduction,  the  latent  heat  liberated  per  unit  volume  of  sample  must  exceed  the 
sensible  heat  required  to  raise  unit  volume  of  the  sample  from  the  undercooling  temperature  to  the 
equilibrium  freezing  temperature: 

wL  -■ISR  .  w  S,)  \()F  0S\  (1) 

where  w  is  mass  of  water  per  unit  mass  of  rock,  L  is  latent  heat  of  fusion  of  ice,  S#  and  S,  are 
specific  heats  of  rock  and  ice  respectively,  and  Op  and  Os  are  equilibrium  freezing  temperature 
and  undercooling  temperature  respectively.  The  test  of  an  observation  is  that 


where  Op  is  the  apparent  equilibrium  freezing  temperature.  Taking  SR  0.2  cal/g  -  °C, 
Sj  =  0.5  cal/g  -  °C,  L  =  80  cal/g,  this  becomes: 


/2.5  x  10' 

\ 


6.25  x  10 


■) 


Op  -  os 


<1. 


(3) 


By  this  criterion,  only  three  of  the  observed  values  for  freezing  point  depression  are  deficient; 
they  are  tests  57,  58,  59  for  saturated  Barre  granite.  However,  it  seems  clear  that  neither  dif¬ 
ferential  thermal  analysis  nor  calorimetry  are  suitable  techniques  for  determining  pore  water 
freezing  characteristics  in  rocks  of  very  low  porosity. 

The  apparatus  used  was  insufficiently  sensitive  to  detect  possible  secondary  phase  changes 
in  the  temperature  range  -35°  to  -40°C,  such  as  have  been  found  in  clays  by  Anderson  and  Tice 
(personal  communication).  The  low  porosity  and  low  surface  area  of  typical  rock  lead  to  water 
contents  per  unit  mass  which  are  smaller  than  those  typical  for  fine-grained  soils,  and  this  means 
that  equipment  for  work  on  rocks  must  be  correspondingly  more  sensitive  than  that  used  for  soil 
studies. 

If  new  apparatus  were  to  be  built,  consideration  would  be  given  to  the  use  of  "long”  cylin¬ 
drical  samples.  This  would  permit  the  results  to  be  analyzed  so  as  to  give  the  thermal  diffusivity 
of  the  rock  relative  to  that  of  the  reference  as  a  function  of  temperature  (see,  for  example, 

Jellinek  and  Fok,  1967). 

Since  unfrozen  water  content  is  almost  independent  of  total  water  content  (Anderson,  1967), 
the  curves  giving  freezing  point  depression  as  a  function  of  water  content  (Fig.  3  and  4)  may  also 
be  interpreted  as  curves  showing  unfrozen  water  content  of  the  rock  as  a  function  of  temperature. 
These  curves  are  used  later  to  check  the  feasibility  of  obtaining  such  data  from  moisture  charac¬ 
teristics  measured  at  room  temperature  in  a  pressure  membrane  apparatus,  as  suggested  by  Keune 
and  Hoekstra  (1967). 

The  logarithmic  plots  of  freezing  point  depression  against  water  content  (Fig.  5  and  6)  indi¬ 
cate  that  for  some  practical  purposes  it  may  be  permissible  to  assume  that  freezing  point  depres¬ 
sion  in  high  porosity  rocks  is  inversely  proportional  to  water  content. 


II 


PHASE  COMPOSITION  Ob'  PORE  WATER  IN  COLD  ROCKS 


I  .at cut  heat  data  lor  sandstone  and  limestone  (Kin.  9- 12)  indicate  that  freezing  is  incomplete 
lot  water  contents  below  about  7  ■  10"  g  water  g  rock.  The  results  accord  well  with  the  finding 
of  Dunn  and  Under  (1965b).  whr  also  found  that  the  water  adsorbed  when  rock  is  exposed  to  a 
1  IK)  relative  Immiditv  environment  is  unfreezable  in  the  normal  sense.  No  abrupt  phase  changes 
weie  detected  for  water  contents  below  the  adsorbed  water  content  (measurements  of  adsorbed 
water  content  are  given  in  Appendix  A). 

An  approximate  calibration  for  latent  heat  can  he  obtained  for  Figures  9  and  10  since  the 
amount  of  water  freezing  at  high  water  contents  can  he  closely  estimated. 

Conclusions 

Low  temperature  differential  thermal  analysis  of  rocks  provides  an  acceptable  method  of 
determining  [tore  water  freezing  characteristics,  although  it  is  not  really  suitable  for  work  on 
rocks  of  vcrv  kiw  porosity,  or  for  work  at  very  small  water  contents.  The  required  measurements 
ate  time-consuming,  and  it  would  be  desirable  to  develop  a  simpler  method  for  estimating  pore 
water  freezing  characteristics  so  that  this  information  could  he  obtained  routinely  in  ordinary 
testing  laboratories. 

For  each  rock  type  at  a  given  water  content  there  is  a  characteristic  depression  of  the  freez¬ 
ing  point  for  initial  nudeation  of  pore  water.  In  the  absence  of  gross  chemical  contamination, 
this  depression  of  the  freezing  point  is  probably  controlled  by  the  pore  geometry.  For  a  given 
rock  type  there  is  a  characteristic  relationship  between  fieezing  point  depression  and  water  con¬ 
tent.  At  high  water  contents,  freezing  point  depression  tends  to  a  limiting  small  value,  and  it  be¬ 
comes  indefinitely  large  as  water  content  tends  to  a  limiting  small  value,  which  appears  to  cor¬ 
respond  with  the  maximum  adsorbed  water  content  of  the  rock. 

If  the  amount  of  unfrozen  water  in  a  rock  is  substantially  independent  of  the  total  water  con¬ 
tent  of  the  rock,  then  curves  relating  freezing  point  depression  with  water  content  can  be  inter¬ 
preted  as  relationships  between  unfrozen  wat«'r  content  and  temperature.  This  is  found  to  be  the 
case  for  frozen  soils  (Anderson,  1967).  By  this  token  the  low  water  content  asymptotic  limit  for 
curves  of  freezing  point  depression  can  be  regarded  as  water  which  is  virtually  unfreezable  under 
typical  environmental  low  temperatures.  Both  the  present  writer  and  Dunn  and  Hudec  (1965b)  find 
that  the  amount  of  '•unfreezable”  water  corresponds  to  the  maximum  adsorbed  water  content  of 
the  rock,  and  since  this  can  be  measured  easily  at  room  temperature  (Appendix  A)  it  seems  to  be 
a  valuable  practical  parameter  (see  Appendix  B  for  an  example). 


B.  INDIRECT  DETERMINATION  OF  PORE  WATER  FREEZING  DATA  FOR  ROCKS 

In  studies  described  in  Fait  A,  the  phase  composition  of  pore  water  at  low  temperature  was 
explored  hv  differential  thermal  analysis  (DTA).  Of  particular  interest  were  freezing  point  de¬ 
pression  curves  yielded  by  DTA,  since  these  curves  may  be  interpreted  either  as  freezing  point 
depression  as  a  function  of  water  content  or  as  unfrozen  water  content  in  saturated  rock  as  a  func¬ 
tion  of  temperature.  However,  with  the  existing  DTA  apparatus  these  curves  are  not  easy  to  ob¬ 
tain,  particularly  for  rocks  of  low  porosity. 

For  more  than  30  years  a  relationship  between  the  free  energy,  or  surface  tension,  of  soil 
moisture  and  freezing  |x)int  depression  has  been  recognized  (Schofield  and  DaCosta,  1938),  and 
recently  a  rigorous  treatment  of  the  relevant  thermodynamics  lias  been  given  bv  Low  el  nl.  (1968). 
hi  connection  with  the  work  by  Low  el  a/..  Keune  and  Hoekstra  (1967)  suggested  that,  as  an  al¬ 
ternative  to  difficult  low  temperature  measurements,  the  unfrozen  water  content  of  frozen  soils 
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might  be  calculated  from  standard  soil  moisture  characteristics  measured  at  room  temperature  in  a 
pressure-membrane  apparatus.  They  calculated  unfrozen  water  contents  for  a  range  of  soil  types 
from  moisture  characteristics,  but  no  independent  experimental  verifications  were  given. 

In  view  of  the  difficulty  and  expense  of  DTA  or  calorimetry  for  rocks,  it  was  decided  that  in¬ 
direct  methods  should  be  investigated. 

Air  penetration  tests 

Air  penetration  measurements  were  made  on  rock  specimens  using  a  standard  soils  laboratory 
pressure  membrane  apparatus  (15-bar  ceramic  plate  extractor,  Soilmoisture  Kquipment  Co.).  An 
pressure  was  supplied  to  the  apparatus  through  a  special  manifold  system  fiom  a  compressor  and 
storage  tank  The  maximum  pressure  of  this  apparatus  (15  bar)  was  really  too  small  to  suit  the  ex¬ 
periment,  but  higher  pressure  equipment  was  not  immediately  available. 

The  rock  samples  were  2', -in. -diameter  disks,  Vi  in.  thick.  For  Berea  sandstone  and  Bane 
granite,  9  identical  samples  of  eacli  rock  were  tested,  and  for  Indiana  limestone  12  samples  were 
tested.  The  samples  were  first  oven-dried  at  110°C  and  weighed  on  an  Ainsworth  balance  to 
0.001  g.  They  were  then  saturated  by  evacuation  followed  by  immersion  in  de-aired  water  while 
still  under  vacuum.  After  surface  drying  they  were  reweighed  to  obtain  the  saturation  water  con¬ 
tent.  The  saturated  samples  were  placed  on  the  ceramic  plate  (which  has  previously  been  satura¬ 
ted  under  pressure),  the  pressure  chamber  was  closed,  and  a  fixed  pressure  was  applied  and  main¬ 
tained  for  two  days.  Pressure  was  then  released,  the  samples  were  rapidly  weighed  to  find  the 
water  loss,  and  they  were  returned  to  the  pressure  chamber.  This  procedure  was  repeated  at  in¬ 
creasing  increments  of  pressure. 

There  were  some  teething  troubles  with  the  apparatus,  which  had  not  previously  been  used, 
and  two  ceramic  plates  were  broken  during  test  runs.  Each  run  suffered  one  interruption,  and  the 
test  results  show  discontinuities  corresponding  to  the  interruptions.  This  is  perhaps  caused  by 
hysteresis  effects,  which  are  known  to  exist  for  soils.  However,  these  small  imperfections  in  the 
data  have  no  significant  effect  on  the  final  results  of  the  experiment. 

The  test  results  are  shown  in  Figures  13  and  14,  where  water  content  and  air  penetration 
volume  are  plotted  against  the  pressure  differential.  It  might  be  noted  that  these  results  can  be 
used  to  make  an  estimate  of  pore  size,  assuming  that  the  pores  are  tubular  capillaries.  From 
elementary  capillary  theory  (Adamson,  1960): 

A P  .  ,4) 

r 

where  AP  is  the  pressure  differential,  y  is  surface  tension  for  water/air  at  room  temperature 
(=  72  dynes/cm),  4>  is  the  contact  angle  (zero  for  water),  and  r  is  pore  radius.  Effective  pore 
radius  is  marked  on  Figures  13  and  14  although  a  logarithmic  scale  is  more  suitable  for  this 
purpose.  . 

Mercury  penetration  measurements 

The  pressures  reached  in  the  air  penetration  tests  were  insufficient  to  provide  freezing  point 
data  for  temperatures  below  -  1.5°C.  Although  high  pressure  air  penetration  apparatus  was  not 
available,  high  pressure  mercury  intrusion  measurements  had  previously  been  made  on  the  three 
rocks  to  determine  pore  size  distributions,  and  there  seems  to  be  no  obvious  reason  why  these 
data  cannot  be  used  to  give  high  pressure  penetration  data  for  water/air  when  adjustments  for 
surface  tension  and  contact  angle  are  made. 
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PORE  RADIUS  ,  micron* 


Figure  13.  Air  penetration  data  lor  Berea  sandstone  and  Indiana  limestone. 


PORE  RADIUS  ,  micron* 


Figure  14.  Air  penetration  data  (or  Barre  granite. 


AIR  PENETRATION  VOLUME,  cm3/g 
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Equation  4  gives  the  following  relationship  between  mercury  penetration  pressure  in  vacuum 
and  air  penetration  against  water  for  a  given  pore  size  r: 

l\Pa  yw 

—  =  —  — -  (5) 

APm  Ym 

where  \Pa  and  \Pm  are  corresponding  pressure  differentials  for  penetration  of  air  and  mercury 
respectively  into  pores  of  a  given  size,  yw  is  surface  tension  for  water/ air,  ym  is  surface  tension  for 
mercury/vacuum,  and  (f>m  is  contact  angle  for  mercury.  Taking  the  contact  angle  for  mercury  and 
rock  surfaces  as  130°,  the  surface  tension  of  mercury  in  vacuum  as  474  dynes/cm  at  25’C,  and  the 
surface  tension  of  water/air  as  72  dynes/cm: 

APa  0.236!>APm.  (6) 

Using  this  relationship,  the  mercury  penetration  data  of  Figures  15  and  16,  which  give. pene¬ 
tration  volume  as  a  function  of  pressure,  can  be  interpreted  in  terms  of  air  penetration  against 
water.  Available  mercury  penetration  data  for  Indiana  limestone  have  to  be  discarded,  as  compari¬ 
son  with  water  saturation  and  air  pycnometer  data  clearly  show  that  they  are  in  error  -  indicated 
penetration  volumes  are  too  low  by  a  factor  of  2.2. 

Calculation  of  unfrozen  water  content  and  freezing  point  depression 

The  object  of  the  calculation  is  to  find  a  relationship  between  water  content  and  pressure. 
Keune  and  Hoekstra  derive  a  relationship  between  soil  water  tension  (or  air  penetration  pressure) 

A P  and  partial  molar  free  energy  AFt; 

AFj  =  AP  (7) 

where  is  the  partial  molar  volume  of  soil  water.  But  the  partial  molar  free  energy  of  water/ice 
in  a  frozen  soil  is  a  function  of  temperature: 

AFl  =  RT  In  (p°/p°)  (8) 

where  T  is  absolute  temperature.  R  is  the  gas  constant,  p°  is  the  vapor  pressure  of  ice,  and  p° 
is  the  vapor  pressure_of  a  standard  state,  chosen  as  bulk  supercooled  water  at  the  temperature  of 
the  sample.  Taking  V{  as  18  cmVmole  for  water,  and  taking  values  of  p°  and  p°  as  a  function  of 
temperature  from  Dorsey  (1940).  the  required  relationship  between  freezing  point  depression  0 
and  pressure  differential  AP  is,  in  the  units  used  for  measurement: 

0  =  5.85  x  10'1  AP  (9) 

where  d  is  in  degtees  Celsius  and  AP  is  in  lbf/in.J  Equation  9  thus  converts  the  measured  rela¬ 
tionship  between  water  content  and  AP  to  a  relationship  between  water  content  and  temperature 
(or  breezing  point  depression). 

In  the  case  of  the  mercury  penetration  measurements  and  their  air  penetration  equivalents, 
penetration  volume  is  given  as  a  function  of  pressure.  For  water-filled  samples:  water  con¬ 
tent/unit  mass  =  (total  pore  vol/unit  mass  -  penetrated  vol/unit  mass)  x  (density  of  water). 

The  porosity  values  taken  in  reduction  of  the  mercury  penetration  data  were  those  measured 
in  the  air  comparison  pycnometer,  viz  Berea  sandstone  20.6%,  Barre  granite  1.2%.  The  dry  bulk 
densities  were:  Berea  sandstone  2.119  g/cm\  Barre  granite  2.636  g/cm*. 
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PORE  RADIUS,  microns 

Figure  15.  Mercury  penetration  data  (or  Berea  sandstone. 
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Figure  16.  Mercury  penetration  data  (or  Barre  granite. 
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WATER  CONTENT  ,  g  wot«r/g  rock 

Figure  17.  Freezing  point  depression  as  a  function  of  water  content  for 

Berea  sandstone. 
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WATER  CONTENT  ,  g  water  /g  rock 

Figure  18.  Freezing  point  depression  as  a  function  of  water  content  for 

Indiana  limestone. 

Comparison  of  calculated  and  measured  freezing  characteristics 

The  freezing  characteristics  calculated  from  air  penetration  and  mercury  penetration  measure¬ 
ments  are  compared  with  those  measured  directly  by  DTA  in  Figures  17,  18  and  19. 

Looking  first  at  the  air  penetration  results,  the  calculated  characteristics  show  little  quanti¬ 
tative  agreement  with  the  DTA  curves;  the  calculated  curve  parallels  the  DTA  curve  for  Indiana 
limestone,  but.  not  fa  Berea  sandstone.  Calculated  freezing  point  depressions  are  less  than  the 
measured  values  for  Indiana  limestone  and  Bane  granite,  but  generally  greater  than  the  measured 
values  for  Berea  sandstone. 

The  mercury  penetration  results  for  Berea  sandstone  give  freezing  point  depressions  which 
generally  are  smaller  than  the  measured  values,  but  the  calculated  and  measured  chatacteri^Mcs 
follow  the  same  trend. 
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Figure  19.  Freezing  point  depression  as  a  (unction  o(  water  content  (or 

Bar  re  granite. 

The  adsorbed  water  content,  which  appeared  to  represent  an  asymptotic  limit  for  the  DTA 
curves  (Part  A)  still  fills  that  role  for  the  freezing  characteristics  calculated  from  mercury  pene¬ 
tration  data. 

Before  any  conclusions  can  be  drawn  regarding  the  usefulness  of  calculated  freezing  charac¬ 
teristics,  the  probable  reliability  of  the  DTA  curves  has  to  be  considered.  The  curves  shown 
give  actual  measured  values,  and  in  view  of  the  consistent  discrepancy  between  measured  and 
calculated  curves  for  high  water  contents  it  is  worth  considering  whether  a  zero  error  in  tempera¬ 
ture  could  have  occurred.  The  temperature  sensing  thermocouples  of  the  DTA  equipment  were 
given  “steady-state"  calibrations  at  0°C,  and  a  zero  check  of  the  apparatus  was  made  on  every 
run  immediately  prior  to  the  first  phase  change.  ;  It  is  possible  that  there  may  have  been  some 
error  due  to  conduction  along  the  thermocouple  wire,  which  was  subject  to  steep  temperature  gradi¬ 
ents  during  cooling  runs.  However,  the  limited  DTA  data  give  no  indication  of  a  correlation  be¬ 
tween  freezing  point  depression  and  cooling  rate,  which  might  be  expected  if  conduction  along  the 
thermocouple  wire  is  a  source  of  error.  The  DTA  records  for  warming  runs,  in  contrast  to  those 
for  cooling  runs,  show  no  distinct  temperature  for  pore  water  phase  change.  They  indicate  that 
melting  begins  at  temperatures  well  below  0°C  and  continues  progressively  up  to  0°C.  In  satura¬ 
ted  specimens  of  sandstone  and  limestone  a  significant  amount  of  melting  takes  place  at  0°C.  The 
other  possible  source  of  error  in  the  DTA  curves  concerns  the  water  contents.  The  water  contents 
plotted  are  believed  to  be  accurate  values  of  the  average  water  content  of  the  sample,  but  at  low 
water  contents  there  is  no  absolute  guarantee  that  the  water  was  uniformly  dispersed  throughout 
the  sample,  even  though  appropriate  precautions  were  taken. 

Taking  into  account  the  possibility  of  error  in  the  DTA  results  and  the  apparent  differences 
between  freezing  and  thawing  events,  the  mercury  penetration  data  give  plausible  calculated 
freezing  characteristics,  but  further  verification  by  calorimetry  seems  desirable. 
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Conclusions 

Air  penetration  measurements  in  a  standard  15-bar  pressure  membrane  apparatus  do  not  pro¬ 
vide  a  good  guide  to  the  freezing  characteristics  of  rocks.  The  maximum  pressure  reached  cor¬ 
responds  to  a  temperature  of  only  -  1.5°C  and,  on  the  basis  of  the  present  tests,  the  agreement 
with  DTA  results  is  poor. ;  However,  with  high  pressure  apparatus  and  improved  technique  it  ought 
to  be  possible  to  obtain  better  results. 

The  freezing  characteristics  calculated  from  mercury  penetration  data  are  quite  encouraging, 
although  they  cannot  yet  be  accepted  as  completely  valid  on  the  basis  of  presently  available  data. 
Equipment  and  technique  for  mercury  penetration  measurements  are  well  developed,  and  the  data 
produced  are  of  broad  interest  in  rock  mechanics.  The  measurements  described  here  were  made 
originally  to  obtain  pore  size  distribution  for  general  descriptive  purposes,  so  that  freezing 
characteristics  were  obtained  without  direct  experimental  cost. 

Mercury  porosimeter  measurements  are  offered  as  an  analytical  service  by  specially  equipped 
commercial  laboratories.  This  is  convenient  when  relatively  few  samples  have  to  be  processed, 
but  results  furnished  by  commercial  contract  for  the  present  study  contained  a  gross  error  and 
some  evidence  of  minor  carelessness. 

It  is  concluded  that,  for  a  simple  first  estimate  of  freezing  characteristics  of  rocks,  water 
adsorption  and  mercury  penetration  measurements  may  be  employed. 


C.  ELECTRICAL  CONDUCTIVITY  MEASUREMENTS 

The  foregoing  investigations  (Parts  A  and  B)  provide  detailed  information  on  freezing  char¬ 
acteristics  of  pore  water  from  0°C  down  to  about  -  15°C,  but  the  methods  used  are  unsuitable  for 
investigation  of  the  state  of  the  pore  water  at  very  low  temperatures,  and  an  alternative  technique 
must  be  sought. 

Although  little  is  known  about  the  electrical  properties  of  common  rock-forming  silicates, 
they  are  believed  to  be  good  insulators,  and  the  conductivity  of  typical  silicate  rock  at  normal 
temperatures  is  thought  to  be  controlled  by  the  water  contained  in  the  rock  (Keller,  1966,  1967). 

If  this  is  the  case,  then  electrical  conductivity  ought  to  provide  a  sensitive  means  of  exploring 
the  amount  and  condition  of  unfrozen  water  at  low  temperatures,  since  the  bulk  conductivity  of 
ice  is  low  relative  to  that  of  water.  Experimental  results  show  that  freezing  of  pore  water  does, 
in  fact,  cause  a  sharp  drop  in  the  electrical  conductivity  of  wet  rocks  (Parkhomenko,  1967;  Dumas, 
1962). 

One  objection  to  the  adoption  of  electrical  measurements  as  a  means  of  studying  pore  water 
over  a  range  of  temperatures  is  that  the  dielectric  behavior  of  water  adsorbed  or  dispersed  in 
porous  media  is  still  imperfectly  understood,  so  that  major  research  effort  is  required  if  results 
are  to  be  truly  intelligible. ;  However,  this  objection  can  be  minimized  if  an  electrical  parameter 
is  regarded  as  a  semi-empirical  index  which  is  to  be  measured  as  a  function  of  temperature  in 
order  to  define  broad  trends  and  discontinuities. 

The  following  measurements  should  be  regarded  simply  as  a  first  step  towards  the  investiga¬ 
tion  of  pore  water  at  low  temperatures. ;  Even  with  very  good  apparatus,  accurate  measurement  of 
all  components  of  the  complex  dielectric  constant  for  wet  rock  involves  appreciable  experimental 
problems  (Scott  et  al„  1967),  major  difficulties  being  shielding  and  correction  for  stray  fields 
and  electrode  polarization  (Hoekstra  and  O'Brien,  1969).  Pilot  tests  gave  direct  evidence  of 
these  problems  when  d-c  measurements  and  capacitance  measurements  were  made,  and  for  the 
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General  properties  of  test 

specimens. 
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Berea  sandstone 
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1.25 
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Indiana  limestone 

g.  go 

14.1 

0.654 

0.005-200 

0.2-  1.0 

lie 

0.90 

500  (bubble  diam) 

1.0 

final  tests  only  the  apparent  a-e  conductivity  at  a  frequency  of  10*  Hz  was  measured  as  a  function 
of  temperature.  This  was  convenient,  as  it  permitted  rapid  measurements  to  be  made  on  samples 
whose  temperature  was  changing  continuously. 

Procedures 

Rock  specimens  were  cut  from  clean  quarried  blocks  with  diamond  tools,  using  only  clean 
water  as  cutting  fluid.  Cure  was  taken  to  avoid  contamination  after  machining,  and  clean  cotton 
gloves  were  worn  when  handling  the  samples.  Ice  samples  were  prepared  by  packing  sieved  ice 
grains  into  a  vibrated  mold,  saturating  with  degassed  distilled  water,  and  freezing  (see  Hawkes 
and  Mellor,  in  press).  General  properties  of  the  test  specimens  are  given  in  Table  V. 

Equivalent  conductivity  was  measured  on  a  Wayne  Kerr  impedance  bridge  (transformer  ratio 
arm  bridge,  type  B221A)  with  an  ESI  a-c  generator-detector  (type  861A).  Running  checks  on  fre¬ 
quency  were  made  with  a  frequency  counter  (TS1  Universal  Counter,  Model  361),  Measurements 
were  made  at  only  one  frequency  (10*  Hz),  since  sample  temperature  was  changing  continuously, 
and  three  samples  were  tested  simultaneously.  In  the  early  pilot  tests  both  conductance  and 
capacitance  were  recorded,  but  it  appeared  that  the  capacitance  readings  were  of  doubtful  value 
for  the  sample  geometry  and  the  type  of  cell  used. 

The  test  specimens  were  2.57  cm  in  diameter  and  5.11  cm  long.  The  ends  of  the  rock  cylin¬ 
ders  were  carefully  primed  with  silver  paint,  which  was  applied  when  the  rock  temperature  was 
about  30  C  to  evaporate  the  butyl  acetate  carrier  before  it  could  penetrate  the  rock  appreciably. 
Aluminum  disks,  which  acted  as  electrodes  and  sample  supports,  were  tacked  on  to  the  primed 
ends  with  conducting  silver  epoxy  adhesive. 

The  Lucite  cell  which  carried  the  test  samples  is  shown  in  Figure  20.  The  electrodes  were 
spring  loaded  to  assure  positive  contact  during  thermal  straining,  and  a  space  for  desiccant  was 
provided.  Each  sample  was  surrounded  by  an  air  space  and  temperature  was  measured  in  this 
annulus  between  the  sample  and  the  cell  wall  (thermocouples  cemented  directly  into  the  samples 
proved  troublesome).  The  cell  itself  was  placed  inside  a  large  Dewar  flask,  where  it  was  supported 
against  the  walls  by  soft  rubber  pads  (Fig.  21). 

For  tests  on  ice,  the  electrode  consisted  of  a  disk  of  0.4  mm  brass  mesh  lightly  frozen  onto 
the  ice,  a  disk-shaped  pad  of  find  steel  wool,  and  an  aluminum  disk.  This  sandwich  was  gently, 
but  firmly,  compressed  by  the  electrode  spring. 

Tests  on  rock  were  made  both  with  the  material  completely  saturated  and  as  dry  as  possible. 
The  rocks  were  saturated  by  placing  dried  samples  under  vacuum  for  24  hours,  then  running  in 
degassed  distilled  water  and  maintaining  the  vacuum  for  a  further  24  hours.  To  prepare  the  rocks 
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Figure  20.  Test  cell  used  for  low  temperature  Figure  21.  Arrangement  of  test  cell  in  Dewar  flask, 
electrical  conductivity  measurements. 


which  here  are  called  “dry,”  air-dry  samples  were  placed  in  an  oven  at  SO'C  for  24  hours,  then 
transferred  to  a  200°C  oven  for  20  minutes,  and  finally  loaded  into  the  cell  with  a  charge  of  anhyd¬ 
rous  calcium  sulfate  desiccant  while  still  hot  from  the  200°C  oven.  This  procedure  does  not  re¬ 
move  adsorbed  water  completely,  but  it  was  the  best  that  could  be  done  under  the  circumstances. 

In  the  early  tests,  before  desiccants  were  used,  it  was  found  that  after  simple  oven  drying  con¬ 
ductivity  fluctuated  considerably  with  changes  in  the  relative  humidity  of  the  electrical  laboratory, 
especially  in  the  case  of  the  limestone,  which  appeared  to  be  hygroscopic. 

The  test  cell  was  cooled  by  nitrogen  evaporating  from  a  pool  of  liquid  nitrogen  (Fig.  21).  The 
cooling  rate  was  adjusted  crudely  by  varying  the  level  of  the  liquid  nitrogen  pool.  The  cell  was 
warmed  by  conduction  through  the  Dewar  enclosure,  and  warming  rate  was  moderated  in  the  low 
temperature  range  by  a  shallow  pool  of  liquid  nitrogen.  Cooling  and  warming  rates  were  not  per¬ 
mitted  to  exceed  20°C/hr.  Check  readings  were  made  with  steady  temperature  after  the  Dewar 
had  been  stored  overnight  in  various  cold  rooms.  The  cell  was  recycled  through  some  parts  of  the 
temperature  range  to  provide  a  check.  At  the  lowest  temperatures  it  was  necessary  to  saturate 
a  covering  pad  of  loose  glass  fiber  (Fig.  21)  with  liquid  nitrogen  in  order  to  cool  the  samples  to 
the  limit. 

Sample  temperature  was  measured  by  copper-constantan  thermocouples,  which  were  encased 
in  beads  of  cement  to  produce  a  thermal  lag  similar  to  the  lag  in  the  interior  of  the  rock  samples. 
Temperatures  were  recorded  on  a  potentiometric-input  strip  chart  recorder  (Moseley  model  7100B). 
Calibrations  were  made  in  an  ice  bath  and  in  liquid  nitrogen,  and  against  mercury-in-glass 
thermometers  at  -10  and  -35°C. 
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Results 

Equivalent  conductivity,  or  dielectric  loss,  is  given  as  a  function  of  temperature  for  Bane 
granite,  Indiana  limestone  and  Berea  sandstone  in  Figures  22,  23  and  24,  respectively.  Data  for 
two  samples  of  fine-grained  polycrystalline  ice  are  shown  in  Figures  25  and  26.  In  the  Berea 
sandstone  and  the  polycrystalline  ice,  conductivity  became  too  small  to  measure  for  the  size  of 
sample  used  when  temperature  dropped  below  about  -125  C.  The  same  data  are  shown  with 
logarithm  of  conductivity  plotted  against  the  reciprocal  of  absolute  temperature  in  Figures  27-30. 

The  bridge  system  was  stable  and  reliable  and  conductivity  readings  were  precise  and  re¬ 
producible.  The  weakest  measurement  was  temperature.  Errors  from  thermal  lag  between  sensor 
and  sample  were  minimized  by  making  warming  and  cooling  runs,  and  by  keeping  rates  of  tempera¬ 
ture  change  low. 

Discussion  of  results 

At  temperatures  above  0  C  the  conductivity  of  saturated  specimens  is  relatively  high:  at 
25  C  the  mean  values  were  4.2  *  10’\  3.5  y  10'!  and  2.8  *  10'*  (ohm-cni)"1  for  Berea  sandstone, 
Indiana  limestone  and  Bane  granite  respectively.  Conductivity  of  the  unfrozen  saturated  speci¬ 
mens  varied  with  temperature  at  rates  comparable  with  typical  rates  for  the  temperature  depend¬ 
ence  of  electrolytes  (Keller,  1967).  The  conductivity  of  saturated  rock  is  commonly  correlated 
with  porosity  by  an  empirical  relation,  known  as  Archie’s  law,*  in  which  conductivity  is  propor¬ 
tional  to  porosity  raised  to  a  power  p,  where  p  lies  between  about  1.3  and  2.2  (Keller,  1966,  1967; 
Parkhomenko,  1967).  Keller  (1967)  gives  representative  values  for  the  empirical  constants  of  the 
modified  form  of  Archie’s  law  for  several  rock  types,  but  when  the  present  data  are  substituted 
into  the  appropriate  expiesslons  the  calculated  values  lor  pore  water  conductivity  areb.fe  »  10'’, 

1  -  10*1  and  10  ■  10'5  (olim-cm)'1  for  Berea  sandstone,  Indiana  limestone  and  Barre  granite  respec- 
i  ivety,  Tin*  itprwud  ut  value*  "rnggi  wi  ilutt  for  rock*  which  tuivr  fliudmi  aniJ  lutunutV  with 
distilled  water  the  specific  surface  area  may  have  to  be  considered  as  a  significant  variable  in 
addition  to  porosity. 

At  temperatures  above  0nC  the  conductivity  for  dry  specimens  is  much  lower  than  that  for 
saturated  specimens,  liy  almost  5  orders  of  magnitude  for  sandstone  and  limestone  and  by  more 
than  3  orders  of  magnitude  for  granite.  If  the  specimens  were  dried  more  thoroughly  using 
vacuum  methods,  conductivity  would  undoubtedly  drop  still  further. 

As  the  temperature  is  lowered  below  T’C,  most  of  the  pore  water  in  saturated  rock  freezes, 
causing  a  sharp  drop  in  conductivity;  conductivity  drops  bv  a  factor  of  approximately  20  in  Berea 
sandstone  and  Indiana  limestone,  and  by  a  factor  of  4  in  Barre  granite.  For  polycrystalline  ice, 
the  conductivity  ot  meliwatei  at  u  C  was  an  oidei  of  iiiagiiitude  highvi  than  the  conductivity  oi 
the  original  ice  at  0  1 C. 

At  approximately  -10  C  the  conductivity  characteristic  for  polycrystalline  ice  (Fig.  25  and 
28?  shows  a  ifiittslltoK  flout  stttfhg  tewpmus**  dtspwtitem*  iuMwfu  U  'and  -IV T  to  a  weak  imm* 
peruture  dependence  between  -10  C  and  -35  C.  The  curve  lor  saturated  Berea  sandstone  (lig. 
24)  shows  a  discontinuity  at  approximately  -8°C,  but  this  break  represents  a  change  to  stronger 
temperature  dependence  for  temperatures  lower  than  -8"C.  A  detailed  plot  of  Figure  22  shows  a 
slight  step  at  -8  C  for  wet  Barre  granite,  but  it  is  doubtful  whether  this  is  significant. 


„l’ s" 
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TEMPERATURE,  °C 


Figure  22.  Conductivity  at  10 4  Hz  as  a  function  of  temperature  for 

Barre  granite. 


TEMPERATURE,  °C 


Figure  23.  Conductivity  at  10 4  Hz  as  a  (unction  of  temperature  for 
Indiana  limestone. 
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Figure  24.  Conductivity  at  10*  Hz  as  a  function  of  temperature 
for  Berea  sandstone. 
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Figure  25.  Conductivity  at  10*  Hz  as  a  function  of  tempera¬ 
ture  for  fine-grained  polycrystalline  ice  (sample  3). 
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Figure  26.  Conductivity  at  JO4  Hz  as  a  function  of  temperature 
for  fine-grained  polycrystalline  ice  (sample  L). 
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Figure  28.  Conductivity  plotted  against  1/ T  (or  Indiana  limestone. 

TEMPERATURE.  °C 


80  0  -80  -120-140160  -180 


1  I  I  1 

\ 


10  ’ 


£  IObl 

E  I 

c. 

O  7 
10 

v 

>  I08  - 

§  1 
Q  a 

z  10 

o 

O  f 

b  10  n  I 


•  1  !  |  !  |  I 

BEREA  SANDSTONE 

frequency  I04  H? 


-200 


Dry 


\Soturoted 

\ 

\ 


10 


L  imit  of 
resolution 


\ 

w 

r 


10“ 
I05 
I02 
I01 
10° 
10  1 


UJ 


10 


I . 

.  1  1  ,  i  i‘  I  1  1  .  1  l 

0  0002  0004  0006  0008  0010  0012  0014 

Vt  (  K 

Figure  2fK  Conductivity  plotted  against  I  T  (or  Berea  sandstone. 


DIELECTRIC  LOSS 


PHASE  COMPOSITION  OF  PORE  WATER  IN  COLD  ROCKS 


29 


TEMPERATURE.  °C 


104 

105 


§  I06 

E 
.c 
o 


O 

o 


10 


o 

3 
Q 


I09 

Iff101 

10'" 

ia12 


40  0  -40 

| — r  !  1  i  7  l 

vxL 

\\ 

J  \60Hr 

Meltwater 


-80 


-120 

T 


l 

!  V. 

Multiple 

points 


-140 

llO4 


POLYCRYSTALLINE  ICE 
■■  Sample  ”J' 

*  Sample  "L" 

frequency  I04H; 


..Coincident 
. .  point 


_1 - 1  i—1  -i _ 1  »  i 


1  .1. 


I03 

I0? 

10' 


CO 

if) 

o 


o 

o : 

U 


10°  w 


UJ 

o 
i  s  - 


I0'1  ; 


10 


10 


,2 


0003  0004  0005 

VT 


jOOG  0007 


(°K  ) 


Figure  30.  Conductivity  plotted  against  1/ T  lor  fine-grained 
polycrystalline  ice. 


At  -35"  to  -40°C  the  ice  characteristic  (Fig.  25  and  26)  shows  adistinct  “knee,”  the  temp¬ 
erature  dependence  of  conductivity  becoming  considerably  stronger  for  temperatures  below  -40'  C. 
There  is  no  direct  reflection  of  this  feature  in  the  curves  for  rocks,  although  there  is  a  converse 
trend  towards  weaker  temperature  dependence  as  temperature  drops  below  -35°C  in  wet  Berea 
sandstone  and  wet  Indiana  limestone  (Fig.  23  and  24).  Curves  for  dry  sandstone  and  limestone 
tend  to  “sag”  to  low  values  at  -35Q  to  -40°C. 

At  -55°C  to  -60°C  the  curves  for  sedimentary  rocks  (Fig.  23  and  24)  show  pronounced 
“knees”  for  both  wet  and  dry  specimens.  There  are  similar,  though  more  subdued,  features  in 
the  curves  for  Barre  granite,  but  they  do  not  show  clearly  in  the  small  scale  plot  of  Figure  22. 

Between  -70°  and  -80°C  there  is  a  small  step-like  discontinuity  in  the  curve  for  ice,  and  a 
most  remarkable  “hump”  in  the  curve  for  saturated  wood  (App.  C,  Fig.  Cl).  There  are  no  corres¬ 
ponding  discontinuities  for  rocks,  but  below  -80°C  the  curves  for  Indiana  limestone  and  Barre 
granite  tend  to  limiting  low  temperature  values  of  conductivity,  with  curves  for  wet  and  drv  speci¬ 
mens  converging  to  the  same  value  of  conductivity  at  -195  C. 

The  curves  for  Berea  sandstone  show  no  evidence  of  a  discontinuity  at  -80  ’C,  but  below 
-105  ’C  the  conductivity  for  saturated  rock  falls  rapidly.  Below  -125  C  the  conductivity  of  both 
wet  and  dry  sandstone  became  too  small  for  measurement  with  the  size  of  specimens  used.  The 
curves  for  polycrystalline  ice  show  a  very  sharp  break  at  -120  C,  with  conductivity  dropping 
abruptly  until  it  became  too  small  for  measurement  at  temperatures  below  -126  C. 
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Recognizing  that  pore  water  freezes  progressively  as  temperature  falls  below  0°C,  Keller 
(1907)  suggested  tliat  freezing  of  pore  water  should  be  regarded  simply  as  an  effective  reduction 
of  the  pore  space  occupied  by  water  for  consideration  of  electrical  conductivity  down  to  -60"C. 
While  this  is  undoubtedly  an  oversimplification  of  the  true  facts,  it  may  give  rough  working  esti¬ 
mates  that  are  adequate  for  some  practical  purposes.  If  effects  of  changes  in  pore  water  salinity 
(which  enters  Archie’s  law  only  as  a  first  power)  are  neglected,  then  the  ratio  of  conductivities 
for  two  different  temperatures  is  equal  to  the  ratio  of  unfrozen  water  contents  for  those  tempera¬ 
tures  raised  to  the  power  p,  where  p  =  2.  By  this  argument  the  drop  in  conductivity  by  a  factor  of 
20  when  pore  water  first  freezes  in  Berea  sandstone  and  Indiana  limestone  would  be  produced  by 
reduction  of  unfrozen  water  content  by  a  factor  of  \/20.  or  about  4.5,  which  is  roughly  the  magnitude 
of  the  reduction  shown  by  DTA  and  penetration  studies  in  Parts  A  and  B.  For  Barre  granite  the 
observed  drop  in  conductivity  by  a  factor  of  4  after  initial  freezing  would  correspond  to  reduction 
of  water  content  by  a  factor  of  2,  which  again  is  not  too  far  from  the  observed  drop  in  unfrozen 
water  content.  At  -60°C  the  ordinate  between  the  extrapolated  line  for  unfrozen  saturated  rock 
and  the  actual  curve  for  the  frozen  rock  gives  factors  for  freeze-effect  conductivity  drop  of  872, 

804  and  78  for  Berea  sandstone,  Indiana  limestone  and  Barre  granite  respectively.  On  Keller’s 
assumption  these  factors  would  correspond  to  reductions  in  unfrozen  water  content  by  factors  of 
29.5,  28.3  and  8.83  respectively  for  sandstone,  limestone  and  granite.  The  actual  unfrozen 
water  contents  at  -60°C  are  not  known,  but  since  adsorbed  water  content  appeared  to  define  an 
asymptotic  limit  for  unfrozen  water  in  the  studies  described  in  Parts  A  and  B,  the  ratio  of  ad¬ 
sorbed  water  to  total  water  provides  a  first  approximation.  The  values  for  this  ratio  are  19.  22, 
and  2.4  for  Berea  sandstone,  Indiana  limestone  and  Barre  granite  respectively. 

Actually  low  temperature  dielectric  behavior  of  water  dispersed  in  porous  media  is  very 
complicated,  being  a  subject  of  high  current  research  interest.  The  present  single  frequency  con¬ 
ductivity  measurements  are  wholly  inadequate  as  a  basis  for  discussion  of  conduction  mechanisms, 
but  they  suffice  to  draw  attention  to  certain  temperatures  and  temperature  ranges  which  may  be  of 
special  interest.  ;  Reviewing  the  significant  temperatures  listed  at  the  beginning  of  this  discus¬ 
sion,  a  number  of  observations  may  be  made: 

0°  to  - 10°C .  Temperature  dependence  for  conductivity  of  impure  polycrystalline  ice  is  high 
in  this  range,  and  certain  mechanical  properties  of  the  same  type  of  ice  similarly  show 
temperature  dependence. ;  Creep  rate  enhancement  for  ice  is  disproportionately  high  (Mellor  and 
Testa,  1969),  and  strength  of  snow,  ice  and  ice-bonded  sand  is  strongly  temperature-dependent 
above  -10°C  (Mellor  and  Smith,  1965).  These  effects  have  been  attributed  to  existence  of  a 
significant  intergranular  liquid-like  layer  above  -10°C,  and  effects  of  typical  contaminants  have 
been  considered  (Mellor  and  Smith,  1965).  The  thickness  of  liquid-like  interfacial  water  in  dozen 
soils  decreases  sharply  as  temperature  falls  from  0°to  -5°C,  and  the  mobility  of  unfrozen  water 
is  apparently  much  reduced  below  -10°C  (Anderson,  1968).  In  rocks  there  is  a  similar  sharp  drop 
in  unfrozen  water  content  possible  below  -10°C  (Parts  A  and  B). 

-35°  to  -40°C.  Low  temperature  differential  thermal  analysis  on  frozen  clay  soils  has  dis¬ 
closed  exotherms  associated  with  phase  changes  in  the  interfacial  water  in  this  temperature 
range  (Anderson,  1968). 

-55°  to  -60°C.  Hoekstra  and  Doyle  (1969)  made  dielectric  measurements  on  frozen  montmoril- 
lonite  at  1010  Hz,  and  found  at  -53°C  a  sharp  break  in  slope  on  the  plot  of  conductivity  against 
reciprocal  of  absolute  temperature. 

-80°C.  Below  -80°C  the  properties  of  interfacial  wafer  tend  to  acquire  some  of  the  charac¬ 
teristics  of  the  solid  state  (Anderson,  1968),  and  the  thickness  of  the  interfacial  film  must  tend 
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to  the  low  temperature  limit  of  about  3  A,  It  is  possible  for  cubic  ice  (Ice  Ic)  to  form  by  conden¬ 
sation  from  the  vapor  below  -80°C  (Kumai,  1968). 

-120°  to  -130°C.  The  glass  transition  temperature  of  water  lies  within,  or  close  to,  this 
range;  the  value  determined  by  Miller  (1969)  is  -111°C  and  the  value  determined  by  Yannas  (1968) 
is  -146°C.  The  mechanical  strength  of  frozen  rocks  appears  to  peak  out  near  -120°C  (Mellor,  in 
press). 


D.  THERMAL  STRAINS  IN  COLD  ROCK 

Following  the  studies  described  in  Parts  A,  B,  and  C,  some  measurements  of  thermal  strain 
were  made  in  the  expectation  that  there  would  be  strain  discontinuities  associated  with  phase 
change  in  the  pore  water.  Initial  measurements  were  made  with  an  optical  extensometer,  but 
it  soon  became  clear  that  this  device  lacked  the  sensitivity  and  stability  needed  for  detection  of 
strain  discontinuities.  A  special  recording  dilatometer  was  then  developed,  and  observations 
were  made  on  three  types  of  rock  with  various  water  contents. 

Measurements  made  with  the  recording  dilatometer  were  rapid  “scanning”  measurements,  in¬ 
tended  mainly  for  the  detection  of  strain  discontinuities.  The  inherent  shortcomings  of  this  rapid 
operating  technique  (chiefly  temperature  gradients  in  the  test  specimen)  were  accepted  for  the 
preliminary  study,  as  it  was  necessary  to  obtain  a  broad  picture  before  proceeding  to  more  de¬ 
tailed  work.  The  results  obtained  prove  the  merits  of  the  apparatus  and  outline  the  general  be¬ 
havior  of  rocks  subject  to  thermal  strain  at  low  temperatures,  but  time  limitations  on  the  overall 
project  restricted  the  scope  of  the  study. 

The  recording  dilatometer  can  be  improved,  principally  by  provision  of  better  temperature 
control.  The  problem  of  low  temperature  thermal  strain  in  rocks  deserves  fuller  investigation,  as 
it  has  an  important  bearing  on  thermal  strains  and  frost  deterioration  in  rocks,  concrete,  and 
ceramics. 

Preliminary  tests 

The  first  measurements  of  thermal  strain  were  made  with  a  Gaertner  optical  extensometer, 
which  consists  of  a  pair  of  10  x  telescopes  (focusing  distance  12  in.)  mounted  on  an  invar  bar  at 
a  horizontal  separation  of  10  in.  The  micrometer  telescope  adjustment  reads  directly  to  0.0001  in. 
The  specimen  was  a  rock  cylinder,  1  in.  in  diameter  by  11  in.  long,  to  which  a  pair  of  photograph¬ 
ically  printed  crosshair  targets  were  cemented  at  points  along  a  generator  10  in.  apart.  The 
specimen  lay  horizontally  in  a  cradle  inside  a  special  temperature  control  cabinet,  the  front  of 
which  was  pierced  with  a  pair  of  “peepholes”  to  permit  sighting  by  the  external  extensometer. 
Specimen  temperature  was  changed  in  steps,  and  each  temperature  was  held  for  3  or  4  hours  be¬ 
fore  a  measurement  was  made.  Results  obtained  for  air-dry  specimens  of  Berea  sandstone,  Indi¬ 
ana  limestone  and  Barre granite  are  given  in  Figures  31-33. 

Although  considerable  previous  effort  by  other  investigators  had  gone  into  development  of 
operating  technique  for  the  optical  extensometer,  the  results  obtained  were  poor,  even  for  obser¬ 
vations  on  metal  bar  (see  also  Berg  and  Maffei,  1968).  It  was  concluded  that  this  optical  extens¬ 
ometer  was  unsuitable  for  the  detection  of  strain  discontinuities,  although  it  was  later  used  to 
provide  an  interim  calibration  for  a  recording  dilatometer. 
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Figure  31.  Optical  extensometer  results  lor  air-dry  Figure  32.  Optical  extensometer  results  for  air-dry 
Berea  sandstone.  Indiana  limestone. 


Figure  33.  Optical  extensometer  results  for  air-dry 
Barre  granite. 


Recording  dilatooeter 

Probably  the  most  satisfactory  kind  of  device  for  thermal  expansion  measurements  at  low 
temperature  is  an  optical  lever  dilatometer,  which  can  be  mounted  in  a  cryostat  and  can  have  a 
sensitivity  up  to  10''  cm  (Bunton  and  Weintroub,  1968).  However,  for  the  present  study  it  seemed 
impractical  to  build  or  acquire  a  suitable  optical  lever  device  in  the  time  available  for  the  work. 
The  requirement  was  for  an  apparatus  that  would  accept  rock  specimens  with  a  minimum  linear 
dimension  of  1  in.  (this  size  being  dictated  by  the  grain  size  of  the  rock),  and  would  vary  temp¬ 
erature  in  the  range  +25°  to  -195°C. 
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The  first  device  to  be  built  consisted  of  a  constant  length  baseplate  and  end  frames,  to  which 
were  attached  a  displacement  transducer  and  spring-loaded  push-rods.  The  specimen  was  a  cylin¬ 
der,  1  in.  in  diameter  by  11  in.  long,  and  it  lay  horizontally  inside  a  heavily  insulated  tubular 
jacket,  which  itself  was  laid  on  rollers.  Vapor  from  a  tank  of  liquid  nitrogen  was  piped  into  tin; 
jacket  to  cool  the  specimen,  and  wanning  was  accomplished  by  heat  conduction  through  the  insu¬ 
lating  jacket.  Temperature  was  measured  by  a  thermocouple  taped  onto  the  specimen,  and  temp¬ 
erature  and  displacement  were  recorded  simultaneously  as  functions  of  time  on  a  3-channel  po- 
tentiometric  recorder.  ;The  first  problems  arose  from  instability  of  the  baseplate  and  end  iVames, 
and  they  were  replaced  by  a  massive  granite  base.  Difficulties  were  then  encountered  in  the  push- 
rod  system  and  the  penetrations  of  the  insulating  jacket,  and  the  system  was  abandoned 

The  final  design  of  the  recording  dilatometer  is  illustrated  in  Figure  34,  and  the  experimental 
assembly  is  shown  in  Figure  35.  The  temperature  enclosure  consists  of  a  metal  Dewar  flask, 

4.5  in.  inside  diameter  by  approximately  15  in.  deep,  to  which  liquid  nitrogen  is  added  as  coolant. 


Adjuttmtnt 

Scrtw  Stolhom  Gag* 


Fund 

Quart!  Rod* 

3  *a.  spactd  120* 


Rock  Sp*cim*n 

4  in.  Ig.  I  in.  dia. 


Figure  34.  Recording  dilutometcr. 
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Figure  35.  Recording  dilatometer.  From  left  to  right:  metal  Dewar  em¬ 
bedded  in  sand  box;  cap  and  plug,  with  suspended  specimen;  2-channel 
recorder;  recorder  power  supply  surmounted  by  thermocouple  ice-bath. 


The  Dewar  is  closed  by  a  heavy  steel  lid  and  a  plug  of  foamed  polystyrene  5.25  in.  deep. 

The  steel  lid  is  surmounted  by  a  platform  which  supports  the  displacement  transducer,  a  Statham 
type  G7A  unbonded  strain  gauge  device.  The  rock  specimen,  1  in.  in  diameter  by  4  in.  long, 
stands  vertically  on  a  stainless  steel  platform,  which  is  suspended  from  the  transducer  platform 
by  three  rods  of  fused  quartz  (4  mm  in  diameter).  A  fourth  (Used  quartz  rod  with  hemispheric  end 
caps  slides  freely  at  the  center  of  the  assembly,  providing  a  push-rod  to  actuate  the  transducer 
as  the  test  specimen  changes  length.  The  length  of  the  suspension  rods  equals  the  length  of  the 
push-rod  plus  the  length  of  the  specimen.  The  device  actually  measures  the  change  in  length  of 
the  rock  specimen  relative  to  the  change  in  length  of  fused  quartz,  which  has  a  coefficient  of 
expansion  (=0.55  x  10‘*  “C'1)  an  order  of  magnitude  smaller  than  coefficients  for  typical  rocks. 

The  parallel  push-rod  and  suspension  rods  automatically  compensate  for  the  effects  of  temperature 
gradient  through  the  plug. 

To  minimize  vibration  or  disturbance  of  the  apparatus,  the  Dewar  was  bedded  in  a  box  of 
Ottawa  sand,  and  the  whole  experimental  assembly  was  set  on  a  rigid  slate-topped  bench  in  a  con¬ 
crete-floored  basement  laboratory. ;  Specimen  temperature  was  sensed  by  a  thermocouple  taped  at 
the  mid-section;  the  specimen  was  wrapped  with  aluminum  foil  to  assure  reasonable  uniformity  of 
surface  temperature  and  minimize  evaporation  or  condensation,  and  it  was  jacketed  with  glass 
fiber  insulation  to  minimize  undesirable  temperature  gradients  and  to  avoid  high-frequency  convec¬ 
tive  temperature  fluctuations. 

Temperature  and  displacement  were  recorded  simultaneously  as  functions  of  time  on  a  2-chan- 
nel  potentiometric  strjp-chart  recorder  (Moseley  7100B),  although  an  XY-plotter  might  be  preferable. 
The  potential  sensitivity  of  the  displacement  transducer  far  exceeded  the  requirements  for  this 
work,  and  the  gain  was  cut  back  to  eliminate  “noise”  in  the  trace  and  to  avoid  frequent  baseline 
shifts  on  the  recorder. 
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Temperature, "f 

Figure  37.  Calibration  data  from  optical  extensometer. 


It  was  intended  that  the  apparatus  should  be  calibrated  by  running  tests  on  a  metal  bar  cali¬ 
brated  at  the  National  Bureau  of  Standards,  but  to  avoid  undue  delay  an  interim  calibration  was 
made  by  running  tests  on  bars  of  aluminum  and  cold-rolled  steel,  which  had  themselves  been 
tested  previously  in  the  optical  extensometer  apparatus  (this  apparatus  was  calibrated  against  a 
brass  rod  supplied  by  the  National  Bureau  of  Standards),  Calibration  data  are  given  in  Figures 
36  -38.  With  a  6-volt  d-c.  supply  to  the  transducer  and  a  recorder  setting  of  10  mv  full  scale  (on  a 
10-in. -wide  chart),  the  calibration  factor  was: 

1  indicated  mv  l.ir>.  1(F4  in.  1.13*  10' 1  mm. 


PHASE  COMPOSITION  OF  PORE  WATER  IN  COLD  ROCKS 


:<6 


Figure  38.  Calibration  data  from  optical  extensometer. 


The  main  fault  of  the  apparatus  in  its  basic  form  was  lack  of  temperature  control. ;  This 
shortcoming  could  easily  be  remedied  for  future  studies.  One  possibility  would  be  to  fit  a  cooling 
coil  to  circulate  nitrogen  vapor  around  the  specimen,  with  an  electrical  resistance  heating  coil 
interwound. 

Test  procedure 

Air-dry  specimens  were  usually  tested  by  loading  them  into  the  apparatus  at  room  temperature, 
adding  liquid  nitrogen  to  the  Dewar,  and  recording  a  cooling  run  followed  by  a  wanning  run. 
Specimens  with  water  contents  close  to  saturation  were  usually  prefrozen  slowly  under  controlled 
conditions  (see  Mellor,  in  press,  for  details  of  freezing  techniques).  They  were  transferred  to  the 
dilatometer  while  still  frozen,  taken  down  to  low  temperature  by  liquid  nitrogen,  and  a  wanning 
run  was  then  recorded.  In  some  cases  the  liquid  nitrogen  was  replenished  before  the  pore  water 
thawed  so  as  to  obtain  additional  runs,  while  in  other  cases  the  liquid  nitrogen  was  replenished 
after  the  pore  water  had  thawed  at  the  end  of  the  first  warming  run. ;  In  the  latter  case,  the  re- 
freezing  of  pore  water  on  the  second  run  occurred  abruptly,  in  contrast  to  the  slow  freezing  which 
occurred  when  specimens  were  prefrozen. ; 

In  most  cases  the  initial  portion  of  the  trace  resulting  from  a  freshly  started  cooling  run  was 
judged  to  be  unreliable,  since  rate  of  temperature  change  was  high  (=2°C/min  at  first)  and  stable 
temperature  gradients  in  the  specimen  and  its  jacket  had  not  been  established.  Rates  of  tempera¬ 
ture  change  also  tended  to  be  high  (=1.2°C/roin)  in  the  mid-section  of  a  warming  run  when  the 
liquid  nitrogen  pool  boiled  completely  away  before  the  specimen  had  reached  a  high  temperature. ; 
Apparent  hysteresis  between  warming  and  cooling  runs  is  believed  to  be  due  to  temperature  gradi¬ 
ents  in  the  specimen;  observation  of  consecutive  warming  and  cooling  runs  ought  to  obviate 
errors  from  this  source.  ; 

During  oooling  runs,  the  rate  of  temperature  change  averaged  1.5°C/min  between  +13°  and 
-13°C,  0.9°C/min  between  -13°  and  -41°C,  and  0.6°C/min  between  -41°  and  -71  C.  At 
the  lowest  temperatures,  rates  were  lower  than  these. ;  During  wanning  runs,  the  rate  of  tempera¬ 
ture  change  averaged  1.2°C/min  between  -71°  and  -41°C,  and  0.7°C/min  between  -41°  and  -13°C. 
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Figure  39.  Thermal  strain  records  for  Berea  sandstone. 


Results 

Results  of  the  tests  are  given  graphically  in  Figures  39-41,  and  a  digest  of  numerical  data 
is  given  in  Table  VI. 

The  first  noticeable  feature  is  that  a  strain  discontinuity,  or  "step,”  exists  for  wet  speci¬ 
mens.  In  Figures  39-41  virtually  all  the  strain  "steps”  shown  are  those  measured  during  warm¬ 
ing  runs.  The  step  data  for  cooling  runs  have  been  omitted  from  the  final  plots,  since  in  most 
cases  the  cooling  rates  were  too  high  and  results  are  confused  by  effects  of  steep  temperature 
gradients  in  the  test  specimen.  However,  the  results  which  were  obtained  leave  little  doubt  that 
the  “freezing  step”  occurs  at  lower  temperature  than  the  "thawing  step,"  and  it  tends  to  be 
steeper  on  the  high  temperature  side. 

In  a  thawing  run  the  main  step  effect  is  spread  over  a  fairly  wide  range  of  temperature:  10  C 
or  more.  In  a  wanning  run,  the  first  evidence  of  melting  effect  s  may  appear  at  a  rather  low  temp¬ 
erature,  say  -20°C,  presumably  because  of  melting  in  the  smallest  pores  which  contain  ice  Pro- 
pessive  melting  effects  are  then  evident  up  to  a  temperature  slightly  below  0°C.  This  tempera¬ 
ture,  at  which  melting  effects  disappear,  seems  to  correspond  to  the  characteristic  freezing  point 
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Figure  40.  Thermal  strain  records  for  Indiana  limestone. 
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Figure  41.  Thermal  strain  records  for  Barre  granite. 

depression  for  the  rock.  The  "strain  reversal  temperature,"  at  which  the  effects  of  melting  just 
balance  the  thermal  expansion  of  the  solid  constituents,  is  approximitely  -7°C  for  sandstone  and 
limestone  at  high  water  contents,  and  approximately  -5°C  for  granite  at  high  water  contents.  ; 

The  present  results  give  only  a  rough  idea  of  the  absolute  magnitude  of  the  freezing  strain. 

For  wet  Berea  sandstone  the  freezing  strain  for  quick-frozen  rock  is  from  about  1  *  10*4  to  6  *  10*\ 
The  corresponding  figures  for  Indiana  limestone  are  about  0.9  *  10'4  to  4  *  10'4,  and  for  Barre 
granite  about  0.7  *  10'4  to  1.3  *  10'\  As  might  be  expected,  there  appears  to  be  an  increase  of 
freezing  strain  with  increasing  water  content. 

The  other  obvious  feature  of  the  results  is  the  dependence  of  coefficient  of  expansion  on  the 
mode  of  freezing  for  wet  rocks  at  tow  temperatures.  In  the  range  of  temperatures  from  about  -15° 
to  -90  C,  the  coefficient  of  expansion  for  specimens  which  were  "snap-frozen”  in  the  dilatometer 
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Table  VI.  Thermal  strain  dal  a  for  cold  rocks. 


Test 

loafer 

content* 

(g  water/g  rock) 

Specimen 
length 
at  20°C 
(in.) 

Coefficient 
of  linear 
expansion 

(°c-‘) 

Temperature 
range  for 
expansion 
coefficient 
("C) 

Strain 
reversal 
temperature 
("C ) 

Remarks* 

Berea  sandstone 

r> 

9  x  i(f2  (S) 

4.055 

11.3  x  10'6 

-25  to  -125 

-7.5 

A(w) 

6 

5  x  10'5  (a.d.) 

4.015 

11 

9  >  1()'2  (s) 

4.032 

13a 

9  x  !«' 2 (S) 

4.001 

11.4  x  10'6 
9.5(1  x  10'* 
7.26  x  io'6 

-20  to  -75 
-75  to  -120 
-120  to  -170 

-4.5 

A(w) 

13b 

9  x  l()‘2  (S) 

4.001 

12.6  x  10‘6 
11.2  x  10'* 

- 15  to  -55 
-55  to  -85 

B(w.c) 

15 

7.3  x  10  J 

4.011 

8.93  x  10'* 
7.97  x  10'* 

+ 10  to  -40 
-40  to  -90 

A(w,c) 

18a 

9  x  10’2(s) 

4.003 

10.05  >  10**‘ 
7.02  x  I0'fc 

-20  to  -90 
-110  to  -170 

A(w) 

ISb 

9  x  lo  2  (s) 

4.003 

12.9  x  10'* 

-15  to  -47 

-6.25 

B(w,c) 

23 

1.22  >  10'2 

4.012 

10.2  x  io‘‘ 
8.21  x  io‘* 

-t  7  to  -30 
-40  to  -93 

(w.e) 

26a 

6.2  >  10'2 

4.014 

9.75  x  10'* 
8.65  x  10'* 

-20  to  -100 
-100  to  -145 

-7.75 

A(w) 

26b 

6.2  x  10'2 

4.014 

11.15  x  10'* 

-15  to  -70 

-4.25 

B(w.c) 

27 

5.35  x  10'2 

4.015 

11.1  x  10'* 

9.66  x  10’* 

-20  to  -70 
-70  to  -120 

-7 

B(w,c) 

30 

7.8  x  10'2 

4.019 

11.7  x  10'* 
9.4  1  x  10'* 

-20  to  -60 
-70  to  -100 

-5 

B(w.c) 

Indiana  limestone 

7 

2.4  x  10‘5  (a.d.) 

4.047 

8 

5.4  x  10'2  (s) 

12 

9.1  x  10'3 

4.032 

16 

6.2  x  10'2 

4.046 

19a 

5.4  x  10'2  (s) 

4.033 

3.76  x  10'* 

-20  to  -40 

A(w) 

19b 

5.4  •  10'2  (s) 

4.033 

6.40  x  10'* 

-20  to  -40 

7.75 

B(w,c) 

go 

9. 1  x  10'3 

4.048 

2.44  x  10'* 

♦  15  to  -85 

(W.C) 

25a 

3.77  x  10'2 

4.049 

2.90  x  10'* 

-10  to  -90 

-6.5 

Atw) 

25b 

3.77  >  10'2 

4.049 

5.21  x  10'* 
4.18  x  10'* 

-10  to  70 
-70  to  -95 

-6.5 

B(w,c) 

28 

4.69  x  10'2 

4.048 

5.79  x  10'* 
4.67  x  10  * 

-20  to  -60 
-60  to  -90 

7.0 

Btw.c) 

♦Symbols: 

A  •  specimen  pre-frozen  slowly 
B  •  specimen  snap-frozen  in  dilatometer 
(w)  •  warming  run 


(c)  ■  cooling  run 

(s)  -  effectively  saturated 

(a.d.)  •  air-dry 
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Table  VI  (Cont’d).  Thermal  strain  data  for  cold  rocks. 


IV. sr 

Water 

(.‘ini  lent* 

(g  wati't  g  rock) 

Specimen 
length 
at  20  C 

(in.) 

Coellicient 

of  linear 

expansion 

(°CM> 

Temperature 
range  for 
expansion 
coefficient 
<°C) 

Strain 

reversal 

temperature 

<°C) 

Remarks 

Bane  panite 

a 

1.8  -  It)"’  (a.d.) 

4.093 

t 

8.15  ■  K)'1  (s) 

4.090 

Ha 

8.11  '  111"'  (s) 

4.089 

3.97  x  10“ 

3.52  x  10“ 

-15  to  -60 
-60  to  -150 

A(w) 

14b 

8.0  -  ID'’  (s) 

4.089 

6.69  x  10“ 

6.15  x  10“ 

-80  to  -50 
-50  to  -80 

B(w,c) 

17 

8.415  v  It)*’ 

4.095 

4.59  x  10“ 

♦  5  to  -114 

(w.c) 

Lila 

8.15  v  to-’  (s) 

4.091 

3.74  x  10“ 

3.23  x  10“ 

-13  to  -120 
-120  to  -170 

A(w) 

80b 

8.6  >  U)“  (s) 

4.091 

5.55  x  10“ 

-15  to  -86 

-6.5 

B(w.c) 

81 

1.8  *  10M  (a.d.) 

4.092 

4.83  x  10“ 

0  to  -1 12 

84  a 

8.6  v  10‘‘  (8) 

4.091 

4.69  x  10“ 

3.81  x  10“ 

-10  to  -95 
-95  to  -155 

A(w) 

84  b 

8.6  *  10*’  (s) 

4.091 

5.98  x  10“ 

-10  to  -80 

-3.0 

B(w,c) 

89 

8.6  x  10‘‘  (8) 

4.043 

5.92  x  10“ 

-20  to  -77 

B(w,c) 

at 

2.6  x  10'*  (s) 

4.090 

5.39  x  10“ 

+  10  to  -97 

B(w,c) 

•Symbols: 

A  •  specimen  pre-frozen  slowly  (c)  •  cooling  run 

B  •  specimen  snap-frozen  in  dilatomeier  (s)  -  effe<  ively  saturated 

(w)  •  warming  run  (a.d.)  -  air-dry 


is  significantly  higher  than  the  coefficient  for  specimens  that  were  slowly  prefrozen,  with  some 
opportunity  for  extrusion  or  redistribution  of  pore  water.  The  difference  for  Berea  sandstone  is 
about  10%  to  30%,  for  Indiana  limestone  about  70%,  and  for  Barre  granite  about  30%  to  70%.  ■ 

The  strain-temperature  traces  are  nonlinear  over  a  sufficiently  wide  range  of  temperature,  with 
slope,  or  expansion  coefficients,  decreasing  as  temperature  decreases. 

Dlaousaion 

Comparing  the  dilatometer  results  with  results  of  differential  thermal  analysis  and  electrical 
conductivity  tests,  there  can  be  no  reasonable  doubt  that  the  discontinuity  in  thermal  strain  is 
caused  mainly  by  phase  change  in  the  pore  water  of  the  rock. 

Similar  freezing  strains  have  been  observed  in  water- saturated  porous  glass  by  Antoniou  (1964), 
and  in  porous  glass  with  lower  water  contents  by  Hodgson  and  McIntosh  (1960).  In  saturated  parous 
glass  the  strain  reversal  temperature  forthawing  of  adsorbed  water  was  -12° to  -16°C,  and  the  strain 
discontinuity  for  freezing  occurred  at  a  lower  temperature  than  that  for  thawing.  There  was,  however, 
a  strain  discontinuity  between  0°  and  -10°C  that  was  not  associated  with  heat  effects.  The  porous 
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Glass  used  in  Antoniou’s  experiments  had  saturation  water  content  of  0.21  to  0.24  g  water/g  rock 
and  B.E.T.  surface  area  of  82.5  to  94.5  mVg;  the  calculated  pore  diameter  was  about  6  x  10*‘  mm, 
and  all  the  water  was  adsorbed. 

The  dilatometer  results  obtained  so  far  have  given  no  indication  of  secondary  phase  changes 
at  low  temperature  which  might  correspond  with  those  observed  in  frozen  clay  at  about  -40°C 
(Anderson,  1967). 

The  observed  freezing  strains  are  a  good  deal  smaller  than  the  theoretical  maximum  freezing 
strains  for  a  water- saturated  closed  system,  which  can  be  calculated  as  follows. 

For  complete  freezing  of  pore  water  in  a  saturated  closed  system  of  porosity  n  and  zero  elastic 
modulus,  the  potential  volumetric  freezing  strain  fpv  is: 


'PV 


vol  ice  -  vol  water 
specimen  vol 


O.Odn 


(10) 


and,  assuming  isotropy,  the  potential  linear  freezing  strain  fPL  is: 


‘PL 


0.09n 

3 


0.03  n. 


(11) 


(Specific  volume  of  ice  is  taken  as  1.09  m’/Mg.)  However,  as  the  rock  is  elastic  there  will  be  some 
elastic  restraint  which  reduces  the  potential  freezing  strain.  The  net  volumetric  strain  (NV  will 
be: 


‘nv  ~  ‘PV  ~  n‘EI  (*2) 

where  fEI  i  .3  the  elastic  volumetric  strain  of  the  ice.  Assuming  that  the  ice  exerts  an  effective 
pore  preso'iro  -m  ihe  rock,  the  pressure  on  the  ice  equals  the  tensile  bulk  stress  in  the  rock 
matrix: 


a 

‘El  =  jjj 


(13) 


where  K ,  is  bulk  modulus  of  ice  in  compression  and  K  is  bulk  modulus  of  the  rock  in  tension. 
Thus:  R 


'NV  1  ♦  n  (Kj/Kr) 
and  the  net  linear  freezing  strain  <NL  is: 


NL  1  4  n(K,/KR) 

Taking  the  compressive  Young's  modulus  of  ice  as  1.27  ■  10*  lbf  in.1  (see  Hawkes  and  Mellor, 
in  press)  and  Poisson’s  ratio  for  ice  0.4  (Gold,  1968),  then  K,  2.12  •  10*  lbf /in.1  For  Berea 
sandstone,  Indiana  limestone,  and  Barre  granite,  values  of  Young’s  modulus  in  tension  for 
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appropriate  ninths  of  strain  art',  respectively,  0.278  ■  10* ,  3.07  x  10",  and  4.56  x  10h  lbf/in,' 
(flaiiepv,  Hawkes  and  Mellor,  mi jNibli shed).  Corresponding  values  of  Poisson’s  ratio  from  the  same 
source  are.  respectively,  0.12,  0.2S,  and  0.34.  The  resulting  values  of  K R  for  Berea  sandstone, 
Indiana  limestone,  and  Bar  re  granite  are,  respectively,  1.33  •  10’,  2.33  •  10%  and  4.75  ■  10'’ 
ltd  in.'  Substitution  these  values  in  eq  14  and  15,  it  turns  out  that  elastic  restraint  has  little  direct 
influence  on  freezing  strain  (Table  V II ). 


Table  VII.  Linear  freezing  strains. 


Fork  type 

Observed 

freezing 

strain 

■  10 ‘4 

Theoretical 

freezing  strains 

'PL  'NL 

x  /O'"  >  10M 

Tensile  failure 

strain  for  wet 

rock  (from 

uniaxial  tensile 

tests) 

x  10'" 

Berea  sandstone 

1  to  t) 

54.7 

54. 1 

8. 1 

Indiana  limestone 

0.0  to  4 

39.0 

33.8 

3.5 

Barre  granite 

0.7  to  1.3 

1.9! 

1.88 

5.4 

The  potential  volumetric  freezing  strain  for  an  unsaturated  rock,  rpv,  will  be  less  than  rpv 
if  the  pore  water  can  distribute  itself  through  the  pore  system  as  freezing  takes  place.  Assuming 
that  pore  space  can  be  fully  utilized. 


fPV 


vol  ice  vo l  pores 
specimen  vol 


-  1.09K%  -  " 


(16) 


where  W  is  water  content  (mass  of  water  per  unit  mass  of  rock)  and  yR  is  dry  bulk  density  of  thu 
rock  in  g  'cm'  or  Mg  in’.  Ignoring  the  rather  small  effect  of  elastic  restraint,  the  potential  linear 
freezing  strain  for  an  unsaturated  specimen  rpL  is: 

«j,L  1  (1.091V >R  -  n).  (17) 


With  this  simple  interpretation  the  freezing  strain  would  reduce  to  zero  at  a  water  content  where: 


W. 


n 

1.09  7r  ’ 


(18) 


Values  of  for  Berea  sandstone,  Indiana  limestone,  and  Barre  granite  are,  respectively,  8.54  x 
10"1,  5.65  <  10"%  and  2.41  x  10’’  g  water/g  rock.  This  effect  of  unsaturation  is  clearly  a  powerful 
one  in  reducing  freezing  strain,  but  the  calculated  values  of  are  higher  than  some  of  the  water 
contents  at  which  substantial  freezing  strains  were  observed.  Thus  it  seems  that  the  pore  space 
cannot  be  fully  utilized,  even  when  freezing  occurs  slowly. 

The  failure  strains  for  these  rocks  in  uniaxial  tension  are  known  from  other  studies  (Gariepy, 
Hawkes  and  Mellor,  unpublished),  and  they  are  given  in  Table  VII.  It  can  be  seen  that  the  observed 
freezing  strains  come  close  to  the  strains  for  complete  tensile  failure  of  the  rocks,  while  the 
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potential  freezing  strains  lot  completely  saturated  rock  frozen  rapidly  far  exceed  the  tensile  failure 
strain  for  sandstone  and  limestone.  It  should  also  be  noted  that  internal  cracking  of  these  rocks 
almost  certainly  occurs  at  tensile  strains  well  below  the  ultimate  tensile  failure  strain. 

In  Figure  42  expansion  coefficient  measured  within  the  temperature  range  15  to  90  C  is 
plotted  against  water  content  for  Berea  sandstone  and  Indiana  limestone.  While  these  data  are 
rather  sketchy,  they  do  show  certain  trends:  coefficient  of  expansion  increases  with  increasing 
water  content,  and  for  any  given  water  content  it  tends  to  be  higher  when  the  rock  has  been  frozen 
rapidly.  These  effects  are  tentatively  attributed  to  partial  relaxation  of  freezing  strain,  since  the 
coefficient  of  expansion  for  ice,  about  50  •»  10''  nC'*  (Butkovich,  1957),  is  much  higher  than  the 
coefficient  of  expansion  for  dry  rock. 

Conclusions 

The  dilatometer  used  for  this  work  is  basically  a  satisfactory  instrument,  although  an  improved 
temperature  control  is  required.  Improved  equipment  would  permit  more  quantitative  data  analysis. 

The  rocks  tested  showed  significant  freezing  strains  when  water  content  was  above  50"u  of  the 
saturation  value.  Freezing  strain  increased  with  increasing  water  content  and  was  highest  for  rapid 
rates  of  freezing.  Observed  freezing  strains  for  specimens  witli  high  water  contents  were  lower  than 
theoretical  freezing  strains  for  saturated  closed  systems,  but  tended  to  be  higher  than  freezing  strains 
calculated  on  the  assumption  that  pore  space  in  an  unsaturated  specimen  could  bo  fully  utilized. 

Observed  freezing  strains  approached  the  tensile  failure  strain,  and  they  were  of  a  magnitude 
which  might  be  associated  with  onset  of  internal  cracking. 


E.  ISOTHERMAL  COMPRESSIBILITY  OF  COLD  ROCKS 

The  phase  diagram  for  ice  (Fig.  -FI)  shows  that  pressure-induced  phase  translormations  can 
occur  when  ice  is  compressed  at  constant  temperature.  If  the  pores  of  a  rock  are  filled  with  ice, 
then  it  is  likely  that  the  compressibility  curve  for  that  rock  will  exhibit  steplike  discontinuities 
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Figure  43.  Phase  diagram  (or  ice  in  the  pressure-temper  a- 

ture  plane.  (After  Whalley,  1969.) 

when  the  frozen  pore  water  is  transformed  from  one  type  of  ice  to  another  type  of  different  density. 
To  check  this  point  experimentally,  isothermal  hydrostatic  compression  tests  were  run  on  specimens 
of  rock  at  a  temperature  of  10  °C.  Measurements  were  made  in  collaboration  with  Dr.  Pieter 
Hoekstra  of  USA  CRREL,  who  was  engaged  in  a  program  of  compressibility  measurements  for  soils. 

Test  method 

The  test  equipment  (Fig.  44)  consisted  of  a  high  pressure  die,  of  the  Bridgman  type,  with  Vt-in.- 
diameter  tungsten  carbide  pistons  driven  by  a  300,000-lbf  Tinius-Olsen  universal  testing  machine. 

It  was  basically  almost  identical  to  equipment  described  by  Stephens  (1964),  but  had  temperature 
control  modifications.  The  actual  apparatus  is  described  by  Chamberlain  and  Hoekstra  (1969). 

Rock  specimens  for  these  tests  were  necessarily  small:  diameter  was  0.450  t  0.003  in.  and 
length  was  0.500  *  0.001  in.  The  specimens  were  prepared  by  taking  a  nominal  '/rin.-diameter  core 
from  block  material,  dressing  it  to  finished  diameter  on  a  lathe  with  a  toolpost  grinder,  and  trimming 
the  ends  oil  a  surface  grinder.  Dimensioned  specimens  were  dried  and  conditioned  (see  Mellor,  in 
press),  and  were  finally  inserted  into  0.500-in.-OD  lead  capsules  at  a  temperature  of  -10°C.  The 
open  end  of  each  capsule  was  closed  by  slicing  and  folding  in  the  excess  length,  and  it  was  sealed 
by  pressure  in  the  testing  die. 

Pressure  was  applied  to  the  pistons  at  a  low  rate  (about  40  min  for  a  cycle),  but  while  this  rate 
was  sufficiently  low  for  dissipation  of  heat,  it  may  not  have  allowed  enough  time  for  complete  phase 
changes.  Piston  force  was  measured  by  an  electrical  load  cell  (Stromberg),  and  piston  displacement 
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HIGH  PRESSURE  DIE 


Figure  44.  High  pressure  die  and  loading  machine. 

was  measured  by  a  linear  potentiometer  (Computer  Instruments  Corp.).  Outputs  from  the  load  cell 
and  the  potentiometer  were  fed  to  an  X-Y  recorder  (Moseley)  to  give  a  continuous  trace  of  load 
versus  displacement.  The  loading  apparatus  was  housed  in  a  room  regulated  to  a  temperature  of 
-10°C,  and  the  temperature  of  the  die  was  controlled  precisely  at  -10°C  by  circulating  glycol  from 
a  constant  temperature  bath.  Die  temperature  was  measured  by  a  glass  bead  thermistor  and  a 
Wheatstone  bridge. 

Calibration  runs  were  alternated  between  test  runs,  so  that  each  test  run  was  sandwiched 
between  a  pair  of  calibration  runs.  For  a  calibration  run  a  specimen  of  gold  was  tested.  The 
compressibility  of  gold  is  low  and  accurately  known,  and  therefore  the  measured  piston  displace¬ 
ments  could  easily  be  corrected  for  effects  of  straining  in  the  die  and  pistons. 

Friction  effects  were  adjusted  out  by  making  corrections  from  a  closed  hysteresis  loop  of  pres¬ 
sure  versus  piston  displacement  for  loading  and  unloading.  The  correction  was  taken  as  one  half 


Hi 
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Hu1  pressure  liillereuee  between  loading  and  unloading  after  the  specimen  had  been  cycled  enough 
times  to  produce  a  closed  hysteresis  loop. 


Test  results 

The  first  set  of  data  was  oht allied  for  Berea  sandstone  in  dry  and  water-saturated  conditions. 

1  tie  testing  equipment  was  controlled  t>y  an  experienced  operator,  the  work  was  directly  supervised, 
and  the  results  appear  to  he  good  (Fig.  45,  16).  Tests  were  also  run  on  Indiana  limestone  and 
Barre  granite,  in  both  dry  and  saturated  conditions.  These  latter  tests  were  made  on  a  ‘‘space 
available”  basis  without  dired  supervision,  and  the  equipment  was  controlled  by  a  new  operator 
of  limited  experience.  The  results  are  shown  in  Figures  47  and  48,  but  they  are  not  felt  to  be  very 
reliable,  paiticiilarlv  in  the  case  of  the  specimens  which  were  supposed  to  be  saturated.  Although 

these  latter  . . us  were  prepared  and  frozen  in  a  saturated  condition,  the  results  suggest  that 

some  of  the  pore  water  may  have  been  lost  subsequently.  Figures  47  and  48  are  given  for  purposes 
of  record,  but  the  results  will  not  be  dismissed  further. 

Discussion  of  results 

Figure  45a  gives  results  for  the  first  and  seoond  runs  on  dry  Berea  sandstone.  The  curve  for 
the  first  run  affords  an  opportunity  to  check  on  the  use  of  compressibility  as  a  means  of  measuring 
porosity.  Taking  the  intercept  given  by  the  tangent  to  the  linear  high  pressure  portion  of  the  curve, 
a  value  V  'V{)  P  /V0  0.796  is  obtained,  where  Fs  is  the  volume  of  solid  material  and  is 

initial  bulk  volume.  This  gives  a  porosity  of  0.004,  or  20.4%,  which  agrees  quite  well  with  the 
value  19.8%  obtained  by  water  saturation  and  the  value  20.6  obtained  in  the  air  comparison 
pycnometer. 

Figure  45b  gives  results  for  the  4th  compression  of  saturated  Berea  sandstone,  any  unfilled  pores 
having  been  closed  by  preceding  runs.  From  Figure  45b,  and  from  Figure  46  in  which  the  same  data 
are  plotted  in  the  form  of  -\V/VQ  against  pressure,  it  is  immediately  evident  that  the  compressi¬ 
bility  characteristic  has  discontinuities.  The  first  step  increase  of  volumetric  strain  occurs  at 
approximately  1  kb;  this  can  be  identified  with  the  transition  from  lce-I  to  water,  which  ought  to 
take  place  at  1.11  kb  at  -10 ’C  (Table  VID).  Since  the  ice  represents  about  20%  of  the  total  volume 
of  the  unloaded  sandstone,  the  step  volumetric  strain  of  2.5  x  10'1  is  of  the  expected  magnitude  for 
the  Ice-1  to  water  transition  (10.2  >  10’1  -  see  Table  VIII).  The  second  step  increase  of  volumetric 
strain  initiated  at  about  5.15  kb;  tins  pressure  is  higher  than  the  stable  phase  boundary  for  the 
transition  from  water  to  Ice  IV  V,  but  it  corresponds  closely  to  the  met aslable  boundary  (approxi¬ 
mately  5.2  kb)  shown  by  a  broken  line  in  Figure  43.  The  second  step  reaches  completion  at  about 
5.5  kb,  and  there  is  a  third  small  step  at  6.5  kb,  which  can  probably  be  identified  with  the  transition 
from  Ice  V  to  Ice  VI  at  6.25  kb  (Table  VIII),  The  total  volume  change  forthe  second  and  third  steps 
corresponds  closely  to  20%  of  the  total  volume  change  for  compression  from  water  at  4.42  kb  to 
Ice  VI  at  6.25  kb.  The  transition  from  Ice  VI  to  Ice  VIII  at  20.8  kb  is  not  reflected  by  a  distinct 
step  in  the  compression  curve  for  Berea  sandstone,  although  the  curve  does  show  increased  com¬ 
pressibility  above  22  kb. 

These  results  are  encouraging,  although  it  would  be  desirable  to  have  more  data,  with  tests 
run  more  slowly  in  order  to  allow  ample  time  for  phase  changes  to  occur  at  the  appropriate  pressure. 
Present  indications  are  that  the  compressibility  and  phase  transformation  effects  of  frozen  pore 
water  can  be  accounted  for  from  the  phase  diagram  for  water. ;  It  has  been  shown  that  compressi¬ 
bilities  for  rocks  can  be  predicted  with  adequate  accuracy  from  known  compressibilities  of  the 
constituent  minerals,  using  Voigt-Reuss-Hill  averaging  methods  (Brace,  1965).  Chamberlain  and 
Hoekstra  (1969)  have  shown  that  similar  methods  can  be  applied  to  predict  the  compressibilities 
of  soils,  and  they  have  succeeded  in  accounting  for  the  influence  of  frozen  pore  water,  in  both 
completely  and  incompletely  saturated  soils,  by  drawing  on  data  from  the  phase  diagram  for  water. 
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Figure  45.  Pressure-volume  relationships  lor  hydrostatic  compression  of  dry  sandstone  and  water- 

saturated  sandstone  at  -10"  C. 
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Figure  46.  Volumetric  strain  as  a  (unction  of  hydrostatic  pressure  at  -10°C  (or  dry  sandstone  and 

water-saturated  sandstone. 
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Figure  47.  Pressure-volume  relationships  lor  hydrostatic  compression  of  nominally  dry  and  nomi¬ 
nally  saturated  specimens  ol  Indiana  limestone  at  -10°C. 


a  Nominally  Dry  b  Nominally  Saturated 


Figure  48.  Pressure-volume  relationships  for  hydro¬ 
static  compression  ol  nominally  dry  and  nominally 
saturated  specimens  ol  Barre  granite  at  -10‘  C. 
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Table  VIU.  Compressibility  of  ice  at  -10°C. 

! C.ilrulatt><l  from  itriitginaii's  data  hy  Chamberlain  and  llnnkstra,  I960.) 


/•base 

Pressure 

(kb) 

In*  1 

0.0 

let*  1 

1.1 1 

Wairr 

t.ll 

Water 

4. 12 

Ion  V 

4.412 

lOt!  V 

0.25 

lot!  VI 

6.25 

loo  VI 

120.8 

loo  VIII 

120.8 

toe  VIII 

49. 1 

Specific 

volume 

(cm’/g) 

Relative 

volume 

(V/V0) 

1.0900 

1.000 

1.0664 

0.978 

0.9544 

0.876 

0.8688 

0.796 

0.80112 

0.7116 

0.7907 

0.7126 

0.751212 

0.690 

0.7012 

0.644 

0.645 

0.592 

0.60 

0.55 

Conclusions  and  recommendations 

From  the  limited  data  given  here,  and  from  the  results  of  related  studies  which  have  just 
been  outlined,  there  seems  good  reason  to  believe  that  the  compressibility  characteristics  of  ice- 
bearing  rocks  can  be  predicted  if  the  compressibility  characteristics  and  modal  distributions  of 
the  constituents  are  known.  Nevertheless,  it  seems  desirable  to  gather  additional  experimental 
evidence,  of  higher  quality  and  for  a  range  of  rock  types,  before  the  question  is  considered  closed. 


GENERAL  SUMMARY  OF  RESULTS 

In  typical  porous  rock,  freezing  point  depression  ranges  from  a  fraction  of  a  degree  at  water 
contents  near  saturation  to  indefinitely  high  values  10°C)  at  the  maximum  adsorbed  water  con¬ 
tent  of  the  rock.  This  implies  that  pore  water  freezes  progressively  as  temperature  falls,  leaving 
a  residual  quantity  of  unfrozen  water.  The  small  amount  of  water  which  remains  unfrozen  below 
-10  C  is  electrically  continuous,  which  suggests  that  it  is  also  hydraulically  continuous,  and 
perhaps  mobile  under  potential  gradients  (mechanical,  thermal,  electrical).  Changes  in  electrical 
conductivity  with  temperature  suggest  that  the  thickness,  continuity  and  mobility  of  the  unfrozen 
interfacial  water  film  decrease  with  decreasing  temperature  down  to  a  temperature  of  approximately 
-125  C.  Discontinuities  in  the  relationship  between  electrical  conductivity  and  temperature  indi¬ 
cate  rather  abrupt  changes  in  the  properties  of  the  interfacial  layer  at  certain  temperatures. 

The  pore  water  freezing  characteristics  of  rocks  can  be  determined  by  differential  thermal 
analysis,  but  the  required  measurements  are  time-consuming,  and  the  method  is  inherently  unsuit¬ 
able  for  measurements  on  rocks  of  very  low  porosity.  It  appears  that  a  good  estimate  of  the  limit¬ 
ing  unfrozen  water  content  for  typical  environmental  low  temperatures  can  be  obtained  from  a  simple 
mom  temperature  measurement  of  adsorbed  water  content.  It  also  appears  possible  to  estimate  the 
[tire  water  freezing  characteristics  indirectly  from  relatively  simple  mercury  penetration 
measurements. 

When  pore  water  freezes,  the  resulting  9"o  increase  in  specific  volume  can  produce  significant 
volumetric  strains  in  the  rock.  Bulk  freezing  strain  increases  with  increasing  water  content  and 
increasing  freezing  rate;  at  high  water  contents  ('50”i  saturation)  freezing  strain  may  be  sufficient 
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to  cause  internal  cracking  of  I  he  roek.  A  sensitive  record  inn  dilatoineter  with  good  leinperatnre 
control  provides  a  sensitive  means  of  studying  pore  water  phase  changes. 

When  ice-bearing  rock  is  subjected  to  sufficiently  high  hydrostatic  pressure  at  constant  temp¬ 
erature,  ice  in  the  |x>res  undergoes  phase*  transformations,  changing  from  Ice-1  to  successively 
denser  modifications  as  pressure  is  increased.  The  isothei'.nd  compressibility  character! sties  of 
ice  bearing  rocks  can  probably  he  calculated  to  iui  aceepiable  degree  of  accuracy  if  the  compres¬ 
sibility  characteristics  mid  modal  distributions  of  the  constituents,  including  ice,  are  known. 

In  view  of  tin*  demonstrated  existence  of  unfrozen  mobile  water  in  cold  rocks,  reappraisal  of 
frost -damage  processes  in  rocks  and  synthetic  porous  materials  seems  indicated.  A  logical  pro¬ 
gression  in  such  an  investigation  would  involve  study  of  water  migiation  under  potential  gradients 
and  water /iee  segregation  in  preferred  locations. 
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APPENDIX  A:  WATER  ADSORPTION 

Water  adsorption  measurements  were  made  on  samples  of  the  three  rocks  in  anticipation  of  a 
correlation  between  adsorbed  water  and  unfrozen  water.  Disk-shaped  samples,  2  7,  in.  in  diameter 
by  approximately  '»  in.  thick,  were  oven  dried  for  II  days  at  110‘’C  and  weighed  while  still  hot. 
They  were  then  placed  on  point  supports,  exposed  to  an  atmosphere  at  100%  relative  humidity  for 
•1  days,  and  reweighed. 

After  the  first  adsorption  measurement  the  test  was  repeated  and  it  was  found  that  more 
water  was  adsorbed  on  the  second  test.  As  a  result  of  this  finding,  it  was  decided  that  further 
repetitions  of  the  test  should  he  made.  The  results  obtained  are  tabulated  below,  and  are 
plotted  in  Figure  Al. 

Surface  areas  of  the  rocks  were  measured  by  the  B.E.T.  method  using  nitrogen  adsorption: 

Berea  sandstone  1.25  mVg 

Indiana  limestone  0.(554  mJ/g 

Bane  granite  0.105  mVg 

Thus  the  water  per  unit  surface  taken  up  by  each  rock  when  exposed  to  a  humid  environment  is 
approximately: 

Berea  sandstone  4  v  10‘5  g  water 'm* 

Indiana  limestone  4  *  10*’  g  water /m‘ 

Barre  granite  10  *  10'5  g  water/mJ 

Assuming  the  densitv  of  adsorbed  water  to  be  equal  to  that  of  bulk  water,  and  taking  the 
thickness  of  one  molecular  layer  of  water  as  3  angstroms,  the  two  porous  rocks  take  on  an 
amount  of  water  equivalent  to  a  uniform  layer  13  molecules  thick.  The  corresponding  figure  for 
the  granite  is  33  molecules  thick.  Since  the  adsorbed  layer  is  not  likely  to  be  more  than  about 
10  molecules  thick,  it  seems  that  some  of  the  water  is  taken  on  by  capillary  condensation  in  the 
smallest  pores,  at  least  in  the  granite,  which  has  a  much  smaller  proportion  of  large  pores  than 
the  other  two  rocks.  Data  on  pore  size  distribution  measured  by  mercury  intrusion  are  given 
below. 


Table  AI.  Water  adsorption  in  100%  relative  humidity  environment. 

Adsorbed  water  <g  water  'g  rock) 

Berea  sandstone  Indiana  limestone  Barre  granite 
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APPENDIX  A 


TEST  NO  (Number  of  wet-dry  cycles) 


Figure  Al.  Water  adsorption  for  successive  exposures. 
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APPENDIX  B:  ADSORPTION  AND  ABSORPTION  BY  ROCHESTER  SHALE 

In  the  two  porous  rocks.  Berea  sandstone  and  Indiana  limestone,  the  adsorbed  water  is  only  a 
small  fraction  of  the  total  amount  of  water  which  can  be  held  by  the  rock  (about  5%  to  8%  in  both 
cases).  In  the  Barre  granite,  the  adsorbed  water  is  about  45%  of  the  total  water  which  can  be 
adsorbed. 

Adsorption  and  absorption  measurements  were  made  on  another  type  of  rock,  Rochester  shale, 
which  is  the  material  underlying  the  caprock  (Lockport  dolomite)  at  Niagara  Falls.  The  test 
samples  were  cut  from  a  fresh  core  drilled  by  the  U.S.  Army  Corps  of  Engineers,  Buffalo  District, 
above  the  American  Falls  at  Niagara.  The  results  were: 

Adsorbed  water:  (A)  1.39  *  10'2,  (B)  1.41  *  10‘2  g  water/g  rock 
Absorbed  (total)  water:  (A)  1.55  *  10‘2,  (B)  1.58  *  10'2  g  water/g  rock. 

In  this  rock,  therefore,  about  90%  of  the  total  water  which  can  be  held  is  adsorbed  water.  This 
adsorbed  water  is  expected  to  be  virtually  unfreezable  under  the  winter  temperatures  normally  ex¬ 
perienced  at  Niagara. 

Dunn  and  Hudec  found  a  strong  correlation  between  the  percentage  of  “unfreezable”  water 
and  the  susceptibility  of  a  rock  to  freeze/ thaw  deterioration.  According  to  such  a  criterion, 
Rochester  shale  would  be  a  prime  candidate  for  rapid  deterioration  under  frost  conditions. 
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APPENDIX  C:  LOW  TEMPERATURE  CONDUCTIVITY  OF  SATURATED  WOOD 

During  the  tests  on  ice,  one  chamber  of  the  cell  was  uncommitted,  and  it  was  used  to  make 
some  measurements  on  saturated  wood.  The  sole  motivation  was  idle  curiosity  Ixirn  of  log-splitting 
tit  low  temperatures. 

The  test  sample  was  a  cylindrical  piece  of  hard  maple,  1  in.  in  diameter  by  l7,  in.  long.  It 
was  cut  from  a  fresh  length  of  dowel,  and  the  ends  were  primed  with  silver  paint.  It  was  soaked 
in  distilled  water  for  3  days  before  the  test  began. 

The  results  are  shown  in  F’igure  Cl.  One  noteworthy  feature  of  the  curve  is  the  peak  which 
occurs  at  -80  C.  There  is  also  a  “step”  at  about  -30'  C  (-40  is  considered  to  he  a  good  temp¬ 
erature  for  log  splitting  -  temperatures  above  -30  are  not  as  satisfactory). 


TEMPERATURE.  °C 

Figure  Cl.  Conductivity  at  10 *  Hz  as  a  ( unction  of  temperature 
for  saturated  wood  I  hard  maple  l 
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